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RAXBOD: A FORTRAN PROGRAM FOR INVISCID TRANSONIC 

FLOW OVER AXISYMMETRIC BODIES 


By James D. Keller and Jerry C. South, Jr. 

Langley Research Center 

SUMMARY 

A program called RAXBOD is presented for the analysis of steady, inviscid, 
irrotational , transonic flow over axisymmetric bodies in free air. Instruc- 
tions on program usage and listings of the program and sample cases are given. 

INTRODUCTION 

The program described in this report is for the analysis of steady, 
invi«‘cid, irrotational, transonic flow over axisymmetric bodies in free air. 

It solves the exact equation for the disturbance velocity potential and uses 
the exact surface boundary condition. Most of the background about the equa- 
tions solved and the difference scheme used is given in reference 1. This 
report gives instructions on the use of the computer program and also some 
additional details which were not given in reference 1. 

The next section gives a general description of the problem and the 
method of solution. Then the instructions for using the computer program and 
a description of the inputs and outputs are given. The appendices contain 
additional details about some specific parts of the program as well as 
listings of the program and the sample cases. 
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GENERAL DESCRIPTION 

One of the important cons1:'!?rations when trying to solve the full potential 
equation is the choice of a coordinate system. For complex three-dimensional 
shapes cartesian coordinates may be best; however, for simpler two-dimensional 
or axisymmetric shapes the use of a coordinate transformation such that the 
body lies along a coordinate line can greatly simplify the application of the 
exact boundary condition at the body surface. The program described in this 
paper uses a body-nomal coordinate system for closed bodies. For open bodies 
(i.e. bodies with a sting or simulated wake) it uses a body-normal system on 
the forebody up to the first horizontal tangent and a sheared cylindrical 
coordinate system aft of that point. This coordinate system is suitable for 
closed bodies which are blunt on both ends and convex and smooth over the entire 
body or for open bodies which are blunt-nosed and convex and smooth up to the 
first horizontal tangent. It is possible to treat pointed bodies and bodies with 
slope discontinuities but the coordinate system is not well-suited for them and 
their solution may not be as accurate as the blunt-body solutions. 

A stretching is applied to both the normal and tangential coordinates such 
that the infinite physical space is mapped to a finite computational space. 

Thus, the boundary condition at infinity can be applied directly and there is 
no need for an asymptotic far-field solution. Details about the stretching 
functions are given in appendix A. 

The general method of solution is to replace the governing second-order 
partial differential equation with a system of finite difference equations, 
including Jameson's "rotated" difference scheme (ref. 2) at supersonic points. 
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POOR QUALITY 



The difference equations are solved by a co !jmn relaxation method. In order 
to get both rapid convergence and sufficient resolution, the relaxation is 
generally done on three different grids. The difference equations are first 
solved on a crude grid (about 400 grid points) which yields rapid convergence. 
Interpolation of this solution is used as an initial condition for a refined 

- 

grid. This procedure can be repeated to any desired refinement within com- | 

puter time and storage limitations. ! 

The boundary condition at the body surface is applied through the use 
of dumny points inside the body. Details of this computation are given in 
appendix B. 

PROGRAM USAGE 

The program was written in the FORTRAN programming language for use on 
a COC 6600 computer operating under the NOS 1.0 operating system at Langley 
Research Center. The program is overlaid in order to reduce the computer 
storage required. One of the overlays uses several subroutines from the 
Langley Research Center graphics library to create a plot vector file which 
can then be post- processed in order to obtain plotted output. Some modifi- 
cations to the program might be required in order to obtain plots on a 
different computer system. 
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The input cards for each case are summarized in the following table: 




X0(I), 1=1, IXY 

Y0(I), I = 1, IXY 

DYOXN, DYDXT, YMAX, XREF 

IMAX, JMAX, MIT, MHALF, KLOSE, NPLOT 

RFl, COVERG, QF3 

DNDYO, ALF, OXIDXO, XM, CXM, OXIOXM 
GAM, AMINF 


Format 


8E10.3 

8E10.3 

8E10.3 


8E10.3 

8E10.3 


8E10.3 


The definitions of these input variables are as follows: 

OESC - Description of case. Up to 80 alphanumeric characters. Appears 

on printed and clotted output. 

IXY - The number of coordinate pairs used to describe the body. Pre- 

sently limited to 100. 

XO - Input coordinates in the axial direction - 8 per card. 

YO - Input coordinates in the radial direction - 8 per card. 

DYOXN - Body slope at the nose. If it is infinite (as it is for 

blunt bodies) put in a value greater than 900. 

DYDXT - Body slope at the tail (with proper sign). If it is infinite 
put in a value greater than 900. 

YMAX - Maximum body radius. Used to calculate the reference area in 
computing the drag coefficient. 
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XREF 

IMAX 

JMAX 

MIT 

MHALF 

KLOSE 

NPLOT 

RFl 

COVERG 

QF3 


I 

I 

( 

1 




f 


Body reference length. Used for scaling plots, XREF will scale 
to 5 Inches on plots. 

Number of grid lines in the tangential direction. I = 1 is the 
forward stagnation line. I * IMAX is the rear stagnation line 
for closed bodies and downstream infinity for open bodies. 

For each grid refinement IMAX is increased such that 
IMAXj^Ew * 2 (IMAXqi^jj) -1. The present limit on IMAX is 81. 

Number of grid lines in the noma'i direction. J = 1 corresponds 
to an infinite distance from the body and J * JMAX is on the 
body. The same formula and limit that apply to IMAX also 
apply to JMAX. 

Maximum number of iterations (complete relaxation cycles) allowed 
on the first grid. MIT is doubled for each grid refinement. 

Number of grid refinements to be done. 

Body type. 

= 0 for open body (i.e. one with a sting or wake). 

= 1 for closed body. 

Plot trigger. NPLOT = 1 causes write on disc for input to plot 
routines and calling of plot routines. 

Relaxation factor for subsonic points. Usual value is about 1.4. 
Should be in the range 0<RF1<2. The program automatically 
reduces RFl by 10 percent if: (1) The maximum correction, averaged 
over 10 cycles, is greater than that for the previous 10 
cycles, and (2) the last maximum residual occurred at a sub- 
sonic point. 

Convergence criterion control parameter. Usual value is 1. 

Iterations stop when the maximum residual is less than 
C0VERG/(IMAX-1)2, This criterion is the order of the finite 
difference truncation error for subsonic points. If this 
degree of accuracy is not required, COVERG can be made larger. 

Supersonic damping factor for improving iterative stability 
(at the expense of a slower convergence rate). Usual value 
is 0.1, but many cases with subsonic free streams are 
successful with QF3 = 0. Definitely need some QF3 on fine 
meshes with supersonic free streams. Note that QF3 has no 
effect on the accuracy of the converged solution, only on the 
stability and convergence rate. QF3 is automatically 
increased if: (1) The maximum correction, averaged over 10 
cycles, is greater than that for the previous 10 cycles, and 
(2) the last maximum residual is at a supersonic point. 
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DNDYO 


ALF 


DXIDXO 


XM 


CXM 


i 


t 


Derivative of the normal coordinate stretching function at the 
body, . The value of DNDYO can be determined by 

choosing the desired step size for the first grid next to the 

body, hr) . Then (^) = 

° 'Qt/y ^ q 

where AY = 1/(JMAX-1 ) . See Appendix A, 

Exponent in the normal coordinate stretching function, a. Usual 
value is 1.3. Larger values of ALF move the last finite value 
of n farther away from the body and smaller values move it 
closer. See Appencix A. 

Derivative of the tangential coordinate stretching function at 
the nose, _ q. Since AX = l/(If^X-l) then 

ACq ^ DXIDXO/ ( IMAX - 1), which can be used to determine what 
value of DXIDXO to use. It is usually best to use 
A^qSsAHq. The above relation for AC^ is only approximate 
however, and it might be necessary to adjust DXIDXO to c t 
the desired AC^. See Appendix A. 

Axial location, x^^, (in physical coordinates) of the junction 
(or matching point) between the two tangential stretching 
functions, for open bodies only. See Appendix A. Usual 
value about the same as the body length. 

Value of the computational coordinate, X, at the matching point 
of the two stretching functions (for open bodies only). Since 
X varies from zero to one, CXM is the fraction of the total 
number of grid points which will be in the first stretching 
region (ahead of x^^). Usual value is about 0.75. 
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I ! 

I I 


DXIDXM 


GAM 

AMINF 


Derivative of the tangential stretching function at 


point, 


X=X, 


. (ac). .as DXI0XM/(IMAX-1). 

A, 




for open bodies. See Appendix A, 
Ratio of specific heats. 

Free stream Mach number. 


the matching 
Used only 


1 .) 

2 .) 

3 .) 


4 .) 


The Program Output Is Described Below: 

Listing of body geometry Input. 

Other input values. 

Computed geometric parameters In tangential direction. 

I - Tangential grid index. 

S - Arc length along reference surface. 

X - Axial coordinate. 

Y - Radial coordinate. 

THET - Angle of reference coordinate surface, 0. For closed 
bodies 0 is the same as the body angle, 0g. For 

open bodies 9 = on the forebody and 0 = 0 on the 

afterbody . 

THETB - Body angle, 0g. 

AK - Surface curvature on closed bodies. For open bodies 

AK is the surface curvature on the forebody and 
2 

AK = - on the afterbody. 

HY 

F - Derivative of the tangential stretch function, 
Computed geometric parameters in normal direction. 


AN 


Normal grid index. 
Normal coordinate, n. 
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5.) 


6 .) 

7. ) 

8 . ) 
9.) 


( d Y \ 

dn ) .• 

3 

GH - Stretching function derivative at half intervals. 


GH(J) 


j + 1/2 

Iteration history. 

IT - Iteration number. 

DPMAX - Maximum 4) correction, max 

ij 

I, J location of DPMAX. 


,IT 


.IT - 1 


^ij '^ij 


10, JD - 
RMAX 


Maximum residual, maxiR. 






where is the 


IR, JR - 
ISUB, ISUP 

RAVG ■ 
RFl 
QF3 
NS 

se:c/cy - 


right hand side of the difference equation { with 

2 2 

AX , AY , etc. in denominator). 

I, J location of R*^AX. 

• Indicates if maximum residual occurred at a subsonic 
or supersonic point. 

Average value of the residual. 

Relaxation factor for subsonic points. 

Damping factor for supersonic points. 

Number of supersonic points. 

Time for iteration cycle. 


10 .) 


Time for iterations. 

Tabulated values of C„ and Mach number on the body. 

P 

Drag coefficient by trapezoidal and Simpson integration of the Cp's. 
Rough plot of Cp along the body. This plot is distorted in the 
axial direction because it is for equal spacing in the compu- 
tational space. The asterisks show the level of sonic Cp. 

Mach number chart of the flow field in the computational plane. 
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Numbers printed are the Mach number multiplied by 100. 
from top to bottom. J values from left to right, 
n.) X and y coordinates of the sonic line. 


I values 
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APPENDIX A 

COORDINATE STRETCHING FUNCTIONS 

The normal coordinate stretching function is: 

- AY 

where n is the physical coordinate normal to the booy and Y is the compu- 
tational coordinate which varies from zero at the body to one at infinity. 

The constant A cont'-ols the physical step size at the body, A = (^)y = o ’ 
and for a given value of A, the exponent a controls the size of the last 
finite value of n. Larger values of a move points farther away from the 
body. 

The tangential coordinate stretching is a transformation between the 
physical arc length along the reference surface, C ♦ aod the computational 
coordinate, X, which varies from zero to one. For closed bodies the 
transformation is 

r = + (X - ^) A + B (X - 

where A and B are determined by specifv’inn(^) ^ and requiring that 

C = C at X = 1 . These conditions give A = - max 'dX/x = q 

2 

and B = 4 (:' - A). 

'max 

For open bodies the tangential coordinate stretching is divided into two 
regions with the physical location of the dividing poi''^, being an input 
quantity. Also input is the value of the computational coordinate at the 
dividing point, X^^^. Since the computational coordinate varies from zero to 
one, X|^ is equivalent to the fraction of the coordinates which are upst»-eam 
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of The stretching function for the region from the nose up to is 
^ X + d2 + 33 0 < X < X^ 

In the region from x. to infinity the stretching function is 


5 = ^ 


[X-X„) (l*Xm) 


< X < 1 


The coefficients in these expressions are determined by specifying 


X = 0 


, and 


X = X. 


and tequiring that ^ and be continuous 


at X * X„. These conditions give 
m ’ 


X = 0 


‘1 VdX 


X = X_ 


70 C^ - 22C2 + 2C3 


- 84 C^ + 36C2 - 4C3 


30C^ - 14C2 + 2 C 3 


^m ' ^l^m 


C, = — • 


C 2 = b - a^ 


2X„b 

^ m 

^3 ' roT 
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APPENDIX B 


APPLICATION OF SURFACE BOUNDARY CONDITION 
IN REGION OF SHEARED CYLINDRICAL COORuINATES 


The boundary condition ir the sheared cylindrical coordinates is 



V - y' U = 0 

(Bl) 

where 

U = 1 + 4 .^ - y • 

(B2) 


V = 6 

(B3) 


and yg is the body slope. 

This boundary condition (Bl) can be rearranged to give: 

• 4 ..) (B4) 


get : 


Let W 2 (1 + = DPO 

so that 

g 4 >y = DPO (B5) 




(1 


y +Yb 


Let 


= w. 


7 ^ 


y + y 


and introduce 4 > = g 4 )Y and to 


g4)Y = W2(l + T't'x) 
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T 



First consider "ordinary" dummy points which lie inside the 


above the axis (i .e. n 
in the following sketch: 


y*-AY 


> - 




or Y 


n=-Y, 


< - AY) as 


body and 
shown 


k 



The values of the potential 
computed by first letting <j)y = 


function at ordinary duitmy points are 
^i, JMAX-1 ' ‘^’i , JMAX+1 
2AY 


which can be put into the boundary condition (B5) to give 


‘i’i,JMAX+l ^ ‘f'i JMAX-1 ‘ -|^DP0 (B6) 

This result can be expressed in the more general form (which will be needed 
later) : 


•^i JMAX+1 " '^3 ^^i JMAX-1 '^4 “^i JMAX "'^5 (B7) 

2AY 

by letting w^ = 1, w^ = 0, and Wg = -y 

In cases where the physical location of the dummy point is below 
the axis, the boundary condition is handled differently. Because the flow 
^ield is axi symmetric, the potential at a point below the axis is the same 
as that for a point an equal distance above the axis, as shown in the 
following sketches: 
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T 




I 



t 


X 




Let Y-| be the value of the computational coordinate at the dummy point whose 
potential is desired. A Taylor series expansion for (p at this point (which 
is the same as 4'^ ^ yields: 


: also 


Eliminate (pyy 
‘^i,JMAX+l 


Yj 

^i , JMAX+1 ■ ‘f’i , JMAX + ^1 2 '*’YY 

AY^ 

*^i, JMAX-1 ' JMAX * 2 

from these equations and solve for 

\ ’ 

° ^y2 *1 , JMAX-1 * \ y' * '> 


to get 



Now since 4>y = 
‘^i, JMAX+1 



this can be put into the form 
i, JMAX-1 '^4 '^i.JMAX " '^5 


where 






) 


If Y is the (negative) value of the computacional coordinate that corresponds 

d 

to the location of the axis, then Y.j = AY + 2Y^. 

14 


I 


Y, can be found from the stretching function. The stretching function 

a 

is n * "X “ ^ which can be expanded in a 

series for small Y to give: 

^ = Y . aY^ . y3 . ^ ... 

A reversion of this series gives 

* ' <-!•>’ - sU6£.|Ji2j-2i + ... 

Putting n ~ " yg if'to This gives the value of Y^. 
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APPENDIX C 
PROGRAM LISTING 


JVERLAV(J£PP¥,(I,0) 

»POCR*M RAXBOOS ( INPUT ^OUTPUT, TAPES" I NpiJT,T*PE6sOUTPUT»T*PE«) 

Ik ****«*«*********•***«««*#•••«******••*•••• ***************** 

* PELAXATIPN SOLUTION OP EXACT EQUATION POP OISTURPANCE VELOCITY 

* potential fop axisynhetpic transonic flow 

* (COORDINATE INPUT VERSION) 

* PROGRANMFD 3Y jerry C, south, JR, AND JAMES D, KELLER 

* 


***»•« *i»********«*********i»****»***#***********»****»*»»***********»***» 

CALL 0VERLAYC5HJERPY, 1,0) 

CA-L OVERLAY(5HJEPRY,2,0) 

STOP 

END 

CVERLAY(JEPRv,i,0) 

PROGRAM ONEf' 

L 

c 

c IMPORTANT, WHEN I»DIMENS10N IS changed, id MUST BE SET EQUAL TO 
C M‘w I-DIMENSION, 

A 

i: final VALUE OF !M*x, AFTER ALL GRID-HALVING IS COMPLETED, IS 

C IMAX(F INAL) «(IMAX( INPUT )-l)*(2**MMALF)+l 

C SIMILARLY FOR Jm\x(FINAL), 

C 

c 18” OIMENSIC OF P-ARRAY MUST BE AT LEAST AS BIG AS IMAX(FINAl) 

C 2ND 01M£N:>tCN OF P-ARRAY MUST BE AT LEAST AS BIG AS JMAX (F INAL ) ♦ 1 

C 19 I-AP'^A 3 DIMENSIONED AT LEAST AS BIO AS Imax(FInal) 

C 12 JsA.nAYS dimensioned AT LEAST AS BIG AS JmaX(FInal) 

c XS !N0 YS ARE SONIC PT, COORDS, NO NEED TO CHANCE DIMENSION UNLESS 

c "ORE Than 396 SONIC PTS ARE EXPECTED ( VERY UnlIkElY), SUBROUTINE 

C sONLIN PREVENTS CALCULATION OF MORE THAN 398 SONIC PTS, 

^ IXY IS The number OF INPUT COORDINATES USED TO DESCRIBE TmE BODY, 

t COM'ION BLuKN CONTAINS 9 ARRAYS DIMENSIONED AT LEAST AS BIG AS IXY, 

c PROGRAM riNLl CONTAINS fl ARRAYS AT LEAST AS BIG AS IXY AND 

C 5 ARRAvo ‘T LEAST AS BIG AS IXY+l, 

C 

C «, ^***,*****, ******. A. *»e*«A**«*. ********»**«*. ***•*****.****•**** 

C 


2 

3 

4 
*> 
6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 
16 

19 

20 
21 
?2 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 
S3 

34 

35 

36 

37 

38 

39 

40 
«1 
92 

43 

44 
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c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


COMHON P(81|62) 

common XB( 6 n*VB(en«CP( 81 ) 

COMMON TMET(6n.THETB(6ni8T(6t)*CT(»n»w!(81) »w2(«l)iM3(ei) 

* »M8(8n,M5(6n#x8P(«n#0Po(ftn»F(8n»**<(en»s{en»u8(en#Ph{en 
COMMON *N(Al),6(8t )tCH(ei)*C8(81)»D(81),Xl (6n»X2(8n#M(6n,HR(an 

lfHRP(8n,HRM(8! )*HRMM(81) 

COMMON XS(800; YS(«00) 

COMMON lOf ANMAX,0 nDV0»YMAX,C0»RM80*J8KP,TSP 
COMMON /8L0K1/ XST 
common /BtOKZ/ P1,R*0 

COMMON /0LOK5/ 1M*X, JM*X,C2»RF1 »DPM, IDP*JDP,RPM, IR, JR,NS»GHl02 
1,AO30,0X3O,0XDY,0VS0,DX2,0V2 iKLO8E 
common /BLOKR/ GM80«C0GM],Tur.«iSu,CP0|K$TAH 
COMMON /BL0K5/ JMi,0Y,n,JSUP,J50N,Cr3,I8UB,lSUP,8UM»p 
COMMON /BL0K6/ XO(100),YOUOO)#XOP(JOO)»XOPPnOO)#XOPPP(tOO)» 

* YOPn00)iVOPPnO01 »YOPPP(lOft) ,SOO(iO05 ,1 XY,OYOXN,OyOXT 

COMMON /BLOK?/ 8MAX,S1>XM,XIm, AA,0XIOXO«OXIDXM, A2,A3»X10»X1 1 «CXM, 

* DXtXlOfXREF 
common /BUOXe/ ALF 

COMMON /BtOKR/ N 
DIMEN8ION DESC(6) 

DATA PI/ 3 . iai5R2<>53S8R79/|RAO/57,2R577P513082J/ 


mrnm warning warning WARNING WARNING warning 

OONWT FORGEY TO CHANGE 10 WHEN I-DIMENSION IS CHANGED,,,, 

A************ «!****•********•****«*••***•*•••*******•**•***•***•*•• 

I0B81 

CAU. SECOND (in 
MRnE(6|270) Tl 


FOLLOWING 0 INSTRUCTIONS ESTABLISH TIME TO START CLEANUP OPERA- 
TIONS, when CPU TIME (INCLUDING COMPILE TIM£) COMES withIN TSAF 
SECONDS OF THE TIME LlMlTi Tl.» ITERATION IS STOPPED ANO CLEANUP 
STARTS, JPAR4MS 18 AN LRC SUBROUTINE THAT RETURNS JOB TIME LIMIT 
IN 0(H), second 18 AN LRC SUBROUTINE USED TO MONITOR TMf CURRENT 
TIME, 


ft 


TSApaJO, 

KTLmI 

CALL JPARAMS(D) 

TLROdl) 

REAO(5,2R0) OESC 
IF(E(iF(5)} 20,30 
IF (NFLOT.eo.n RETURN 
STOP 

CONTINUE 
REAO(S.300nXY 
ftEAO(5,350)(XO(l)»l»l,IXV) 
REAO(5,330J(YO(I),lal,IXY) 
REAO<5,330)(^YOXN,DYOXT,VMAX,XREF 
READ(S,500 JIMAX,JMAX,h1T,MHALF,KL08E,NPLOT 


VMIGLNAL page 13 
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as 

48 

<17 

tt6 

<19 

50 

51 

52 

53 

54 

55 

56 

57 
50 

59 

60 
61 
62 

63 

64 

65 

66 
67 
66 

69 

70 

71 

72 

73 

74 

75 

76 

77 
76 
79 
60 

fit 

82 

63 

84 

65 

66 
67 
66 
69 

90 

91 
^2 

93 

94 

95 

96 

97 
96 
99 

100 

lOl 

t02 

103 

104 
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READ(S«330)Rri,COVCRC|OP3 lOS 

REAO(5»33())nNOvO»ALP»OXlOXO«XH»CXM»OXlDXH 106 

RE*D(S»330)GAM, amine 107 

>>RITE(b>2BO}OE8C tOO 

MRITE(6»720)(!,XO(n>rO(I)^!al,lxy) 100 

WRlTE(6f «20)OYDXN,nVOXT#VMAX»XREF 110 

T8AFaTSAF«NPLOT*30, »»l 

THETAN*ATAN(0V0XNJ 112 

TMETATsATAN(oyDXT) 113 

DYDSNaSlN(THETAN) 110 

0X08n«C08(THETAN) 115 

OYOSTaSlNtTHETin ll6 

OXOST«COS(THETAT) 117 

lF(OYOXv.tT,<>00, 100 TO 31 118 

DYOSNSl, 119 

OXDSNaO. 120 

31 CONTPJUE 121 

IF(APS(OYOxn ,UT,900,)G0 TO 32 122 

OYOSTa.l, 123 

0X03Ta0, 120 

32 CONTINUE 125 

CALL FIT ( lxy,xn»YO.SOO,XOP,XOI»P,YOP, YOPP.OYOSN,OXOSN»DYt)STf D*DST ) 126 

CALL SPLlFn»lXY»8O0»X(J|X0?,X0PP,XOPPP,»,OXt)SN»l,l)X08T»INC> 127 

CALL SPLlP(l»lXY|SOO*YO,YOP|YCPP»TOPPP»l*OYDSN,l,OYOSTr)''D) 126 

NHALFaO 129 

ANM*Xal,FY08 1*0 

JSXPal 131 

JPAGEaSl 132 

Nao 133 

00 IF CJMAX/J8KP,LE,JPAGE> 60 TO 50 130 

JSKPaJSKP41 15S 

GO TO i)0 136 

50 CONTINUE 137 

XldJeO, 136 

xamaXHl) 139 

RHS0aYNAX**2 

GMlaGAMal, 101 

GM102a,S*GMl 1«2 

COGMlaGAM/CHl 103 

AMSQ:AMImF«*2 100 

GMS0s6m102*AM3(} 105 

AOSOaGM102tl ./AMSQ 106 

T0GM80a2,/(CAH*AMSO) 107 

PSTARa(2,*U .♦GMSO)/(1,*GAM))«*GOGM1 108 

CP8TARsT0GM30*(PSTAR«1,) 109 

CP0aT0CM9Q*( (1 .♦0MSQ)**G0 Gh1-1 ,) 150 

KSTARa0,S*30,*(CP0-CPSTAR) 151 

IF (KSTAR.GT.IOO) KSTAHelOO 152 

60 CALL SECOND (Tl) 153 

MRITE(6,320) IMAX, Jm^x,miT,MMALF,KLOSE»nplOT 150 

* |RF1 .COVERC, OF3fONDYO,ALF, 0X10X0, XM,CXM,DXI0XM, GAM, XhiNF 155 

156 

****«**t««******<i********«**«*«*«***«****«**A*******A****«****«**« 157 

156 

OVERLAY(Wl) SETS UP THE TANGENTIAL COORDINATES 159 

160 

*«*****«**A«*«*«****i *********** ***•*«*•***«*•*••*• *<»•****•••>«•>••« 161 

162 

CALL 0VERLAY(5HJERRY,1,1,6HRECALL) 163 

C 160 
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A***************************************************************** U5 

U* 

0VEHLAY(1»2) CALLS NTRANF AND W1M2 1*7 

168 

A***************************** ******************* **«•*•••***•*•**• 169 

170 

CALL 0VERLAV(SHjERRY«t*2#6HRfCALL) 171 

0X8O*l,/OX**2 

RCHCRalOO.AOXSQ 

OKOV*,5/(OXaOV) 

0YS0«1,/DY*a2 >75 

0Y2«,5/0V >76 

0X2*, 5/OX >77 

JMlBjMAX-1 >78 

KP0INT*(I*'AX»1)*(JMAX-1) >79 

POINTSoKPOlUT >»0 

WRITC(6,«70) >«> 

00 90 1*1«I»1AX >62 

LS(1)*0 >«3 

TD*THET(I)aR*0 >6R 

TB0«TM€TB(T)*RA0 >65 

WRITEC6, 480) I , S ( I) » X0 ( I ) , Y9 ( I ) , TO, TBO, AK ( I) ,F (I ) 166 

90 CONTINUE >67 

WRITE(6,45i)ALF >6« 

WRITE(6,450) >6’ 

WR1TE(6, 4o0) (J,AN(J),C(J),GH(J),Jsl,JMAX) >90 

CALL secoNO (n >’> 

T»T«T1 >^2 

WRITE(6,430)T >*>3 

IF (NHALF.CT.O) go to 100 >94 

CALL fSTIM (P,I0,INAX,JMAX) >95 

100 IT*0 >96 

00 110 1*1,11 >97 

110 OPOdloSTd) >98 

IF (KLOSE.EO.n CO TO 150 >99 

I2.I1A1 200 

I3*IMAX-1 201 

DO 120 lsI2,I5 202 

120 0P0d)e«(2(l)*(l,+Fd)'‘DX2*{Pd + l,JMAX).pd-l,JMAX))) 203 

0PI*3,*P(IMAX, JMAX)*4,*P(IHAX«1 , JMAX)aP(1MAX-2, JHAX) 204 

0P0(IMAX)*nt2dMAX)*d .♦FdMAX)A0X2*0Pl) 205 

130 00 140 I«1,IMAX 206 

140 P(I, JHAXAn*w3(n*Pd» JHAX-n*K4d)*P(l, JMAX)-w5(I)*0P0d) 207 

MR1TE(6,4'^'0) 208 

CALL SECOND (TO) 209 

8UNlal,e«07 210 

8UH«0, 211 

COVR*COVERfi/FLOAT(IMAX«l)**2 212 

150 CALL SECOND (Tl) 215 

JSUPaO 214 

IF (•MINF.GE.l ,) JSUP*1 215 

JS0N*0 216 

IF (ABS(Af11NF»l,),LE,>.E«'06) JSON*l 217 

218 

•«*«*<.**••*****• A A************************* *•*••«•***«•*•***• •«•*• 219 

220 

0VCRLAY(1,3) is the »«IX£0 FLO* POTENTIAL ITERATION LOOP 221 

222 

223 

224 
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CALL 0VERLAr(5HJCRf)y,l,3,*HRECALLJ 
IT«1T*| 

RAVG*3UMRP/P0INT8 
CALL SECOND (T) 

TUT-Tl 

NRlTe(6»50O)nf0PM#JDP,JOP,RPM,lR,JR, ISU0»I5UP,RAVC»RFUOP3,N8.T1 

c 

c •**«•*«••**«»•** *«**********«««****^«.*«****«. *«•*••*•*«••.•.••• 

c 

c CHICK FOR TJHE LIHIT, 

C 

c •*•**•*****•«****•*•«******«•**•*•*••*••*•*****•.*•**««•*•**•«••* 

C 

IF (TL-T.GT,T8AF) CO TO 160 
HRJTC(6*570)T,TL»PPH»COVR 
KTL«2 
60 TO ISO 

160 CONTiNUt 
C 

C 

J* 

C CHECtf FOR OTVERCfNCf, 

C 

C 

C 

IF(RPH,LT,RCHEK) CO TO ui 
wRITE(6*610) 

GO TO 10 

161 continue 

c 

c «*••«*••*••**•*****••*****•••**«••**•*•«*«***••••**•****«••*•.«•*, 

c 

C CHECK TOR CONVERGENCE OR ITERATION LIHIT 

C 

c 

1F(RPh.GE,C0VR1G0 to 171 
HRIT£(6f700)RPN,C0VR 
CO TO 160 

171 CONTINUE 

IFUT ,LT,HIT)G0 to 172 
hHITE(6,310)HIT,RPm,COVR 
GO TO ISO 

172 continue 
C 

c «*««**«***««*«*•**•*•*•****•***•«•*•*•*«*••***«****•**•**«*.**•••• 
c 

C INCREASE PM1*ST OAnPINC COEFFICIENT OR DECREASE RFl IP AVERAGE 

C maximum correction OF LAST 10 CYCLES MAS INCREASED OVER PREVIOUS 1 

c 

c 

SUMuSUM+OPM 

IF (MODdT.lO) ,NE,0) GO TO 150 
IF (SUMl.CT.SUM) 60 TO 175 
QF3«0F5*. 1*ISUP 
RFl*RFl*(l,«,l*15Ub) 

SUMl«i ,E407 

IFdSUP.Ml.nMRITElb.seOlOFJ 

lF(ISuH,EO,i)f'RITE(6,710)RFl 


225 

226 
227 
226 
229 

250 

251 

252 

255 

256 

255 

256 

257 
256 
259 
2C0 
261 
2C2 

265 

266 

265 

266 
Z67 
266 
269 

250 

251 

252 

253 
250 

255 

256 
25/ 
256 

259 

260 
261 
262 

265 

266 

265 

266 
267 
266 

269 

270 

271 

272 

275 

276 

275 

276 

277 
276 
279 
260 
261 
262 

265 

266 
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60 TO 17« 26S 

17S 8U*«|«8UM 286 

17« SUMaC. 287 

60 TO 150 288 

180 CALL 8ECONO (T) 289 

290 

A***************************************************************** 291 

292 

THE EOLLOaiNO STATEMCNT8 CALL fOH PREPARATION ANO PRINTING OP 295 

CP» MACH NO,, drag, rough CP PLOT, HACH NO, CHART OP FL0*< Pl?LC, 299 

SONIC LINE CALCULATION, AND WRITING ON DISC POR CALtOMP PLOTS, 29S 

296 

******«*•*«*«•**••*****«>•*•****•**•*•******•*****•*****«*••***•*** 297 

298 

TlaT«T0 299 

WRITE<6, 510) T1,IT,NMALP SOO 

501 

**t«********«««*«**«***************o*«**************************** 502 

505 

0VCRLAV(1,<i) calls CP800Y, DRAG, AND CPPLOT, 509 

50S 

****•***<**••*****•*«**••*«•**«** '.••***«.*tt********«***«********** 506 

507 

CALL 0VCRLAV(5HJERRV,l,9,6HRrCALL) 508 

509 

A***************************************************************** 510 

511 

0VERLAV(1,5) CALLS HCHART AND SONLlN, 5l2 

515 

fc*******A«**t«<MtA«*******«ft<t*|i***«*A* <****« O************** «••••*••• 519 

515 

call 0VERLAY(5HJERRY,1,5,6HRECALL) 516 

CALL second (T1) 517 

TlsTl-T 518 

WR1TE(6,600)T1 519 

IF (NPL0T,nE,1) go TO 220 520 

HRITEi4) OESC 521 

wRITE(9) IHax,Jhax,IT,KLOSE,n 522 

WRITE(9)CPSTAR, AM1NF,0PH,XR£P,OXIOXO,0n0Y0,QP5 525 

IP(KLOSE,EO,1)GO TO 211 529 

HRITE(9) CXH,XH,X1M,0X10XH 525 

211 CONTINUE 526 

WR1TE(9HAN(J),J*2,JHAX) 527 

HRITE(«HST(l),Itl,lHAX) 528 

hRITE(9)(CT(1),I81,IhaX) 529 

WHITE(9) (xe(I),Isl,IHAX) 550 

HRITE(O) (YB(I),I«1,1HAX) 551 

NRITE(9) (CP(I),I«1,1hax) 552 

IF (M,e0,0) GO TO 220 555 

NRITE(O) (XS(I),Iai,N) 559 

WRITE(«) (VS(I),1*1,N) 555 

220 CONTINUE 556 

lF(KTL,t0,2'00 TO 20 557 

IF (NHALP,f Q, hhALP) go TO 10 558 

NHALP>NHALP«1 559 

590 

A************** **••**•*•**••*•••****•*•«*«**•»**•••*••*••••*••**•• 591 

592 

0VERLAY(1,6) 18 The grid REPINEWENT routine, 595 

599 
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c ««*««««4**«*««**««*««****««*.«««*******<i**«*ft*«**»**************** US 
c 3«6 

CAU. 0Ve»LAY(5HJE«RY*l»6,6HREC*U) 3«7 

OX*, 5*0* 50^ 

250 If (JMAX/JSKP.lE.JfAGE) GO TO 260 350 

JSKPaJSXf^t 3S1 

GO TO 250 352 

260 COnMNUE 353 

CO ro 60 356 

270 fOPHAT computing TIM£e,f6, J ,0H SECONOS/) 355 

200 format (IHI.BAIO) 356 

200 format (BAtO) 357 

JOO format (t6I5) 350 

JIO F0RMAT(/a •—•010 not converge 1N*I«* cycles, •••• Rmaxba 35R 

• E9,2*, C0VR«*E9,2/) 360 

320 F0RMAT(6HUMAXan/6M JMAK*n/5H MiTslR/TH mmalF»H 361 

• /7H MtOSE*ll/7M NPL0T«II/5 m RH*r5,3 362 

« /OH C0veRGeE9.2/5H QF3«C9,2/7 h 0N0Y0«E10,5 563 

• /SH ALF«F«,2/BH OXlOX08tlO,3/«H XM«E10,3 366 

• /5M CXm«E10,3/6m nxl0XM*El0,3/5H (;amsf«,2 365 

• /7H AminF*F6,0) 366 

330 format (Btl0,3i 367 

370 FORMAT!/* must STOP ITERATIONS, CLOSE To UmE LIMIT,*/ 368 

• • COMPUTING TIME e*F6,l* TIME LIm1Ts*F6,1/ 369 

« • RMAXa*E9,2«, C0VRa*E9,2) 370 

920 FORMAT!/* OYOXN»*F 1 0 , « # /* 0YDXT«*F 1 0 ,9 , /* YMAX»*F 1 0 , «, /• XREFs* 371 

• »F10,«) 372 

930 format !/»9um CPU SECONDS FOR flOOY GEOMETRY CQmPuT A T 1 ONSt , F 6 , 3/ ) 373 

950 format !10X, JHJ,ux,2MAN,10X,lMG,nx,2Mr,M/l 379 

951 FORMAT!/, *l—» NORMAL COORO, S ETCH FOR ALFe*F6,3* •••••*/) 375 

960 format !I12,3E12,9) 376 

970 FORMAT !1M1,9x,1m1,9X,IMS,11X,IHX,11X,1HY,10X,9HTh£T.6X,5mTmETB,B 377 

1X,2MAK,10X,1hF//) 378 

980 FORMAT !112,BE12.9) 379 

990 FORMAT ! IMI ,2X,2H IT, 3X5M0PM AX, SX2mI0,2X2MJ0»3X9HRMAX,6X2h1R,2X2mJR, 360 

1 1X9H1SUB, lX9HlSUP,3X9HHAyG,6X3MRFl ,4X3MOF3,6X2 hnS, 381 

2 3X7MSEC/CYC/) 382 

500 FORMAT!l5,En,3,21«,EU,3,2I9,2I5,Ell,3,2F7,3,l6,F9,S) 3B3 

510 FORMAT!13m0CPU 5EC0n08o,F 7,2,9M F0R,19,19M ITERATIONS, NmaLF»I1/) 389 
600 format 197H0CPU SECONDS TO C-OmPUTE ANO PLOT CP ANQ MCH ARTs , F7 , 3/ ) 385 

610 FORMAT!///* •—•••••••••••••••DIVERGENCE, Rmax EXCEEDS RCmEk,* 3B6 

1 * GO DIRECTLY TO *A1L, 00 NOT PASS 60, DO NOT COLLECT $200, • — * 367 

2 ••••—••••••••••••ft^//) 368 

680 FORMAT!/* OF3 INCREASED T0*F6,3* BECAUSE lO*CTCLf AVERAGE OF* 389 

1 • rmax increased,*/) 390 

700 FORMAT!/* •— CONVERGENCE--^^. RMAX«*E9,2*, C0VHa*t9,2/) 391 

no FORMAT!/* RFl DECREASED T0*F6,3* BECAUSE 10»CYCLE AVG FOR* 392 

1 « RMAX INCREASED,*/) 393 

720 FORMAT!/* INPUT COOROINATE5*/9X1M19XIHX9X1MY/H5,2F10,6)) 399 

END 395 

subroutine F1T!N,x,Y,S,X1,X2,Y1,Y2,OY1,OX1,OY2,OX2) 396 

DIMENSION X!1),Y!1),8(1),K1!1),X2!1),Y1!1),y2!1) 397 

RCS91.0E-07 398 

TOL«,0625*R£8 399 

K«0 900 

KMAXa500 901 

8(1) a 0. 902 

M a N ol 903 

DO 22 lai,M 909 


22 




VAL ■ X(l*n -xd) AOS 

dum ■ Y(i*n •Y(n «06 

22 sn^n • 3(1) ♦SQRT(VAL**2 yOUM**;) «07 

31 CALL SPLl(^(l»^iS,X,Xl,X2,X2,l»DXdl>t>X2,I»JD) 403 

call SPLlPn .N,S,Y»Vl,V2,Y2,l»0nil,0Y2,lND) 409 

ERR ■ 0, 410 

PUW ■ 0( 411 

ro 32 1«1,M 4]2 

SO ■ S(H1) •OUM «}} 

OUM ■ 3(l*n 414 

31 ■ 8(I«n •3(1) 41S 

X3 « (X2(I + 1) •X2(in/Sl 416 

V3 ■ (Y2(Ul) -Y2(m/Sl 417 

CALL ARCL (Sl,SO,Xl(n,X2n),X3,Vl(I),v?n),Y3,R,lNO,TOL) 416 

VAL ■ ABS(S1 -80) 419 

ir (VAL -fRR) 32,32*33 420 

33 ERR « VAL 421 

32 S(I«1) s HI) ySI 422 

KSKtl 423 

1F(K,L£,KMAX)G0 to 34 424 

HRITE(6,9901) 425 

RETURN 4^6 

34 COnTImJE 427 

IF (CHR -RES) 41,41,31 426 

41 RETURN 429 

9901 FORMiT(* FIT FAILEO TO CONVERGE*) 430 

End 431 

SUBROUTINE split («,n,S,F,FP,FPP,FPPP,KM, VM,kn, VN, InO) 452 

C SPLINE FIT - JAMF.SON 453 

DI^'EnSION S(l),F(n,FP(l),FPP(l),FPPPtl) 434 

IND s 0 435 

K • IA9SCN -H) uli, 

IF (K -1) 81,81,1 437 

IK e (N .M)/K 436 

1 ■ R 439 

J C M 4K 440 

OS » S(J) -8(1) 441 

D e OS 442 

IF (OS) 11,61,11 443 

11 OF ■ (F (J) -F(I))/08 444 

IP (KM .2) 12,13,14 445 

12 U ■ ,5 446 

V ■ 3,*(0F -VM)/DS 447 

00 TO 25 446 

1 3 U s 0 , 444 

V • VM 450 

CO TO 25 451 

14 U B • 1 , 452 

V o -DS*VM 453 

CO TO 25 454 

21 : • J 455 

J • J *K 456 

OS s S(J) -S(I) 457 

IF (0*DS) 81,81,23 456 

25 OF « (F(J) -fCin/DS 459 

8 • I, /(OS *03 *U) 460 

U a 8*00 461 

V t B*(6,*0F -V) 462 

25 FP(I) B U 463 

FPPU) B V 464 
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g 

■ 

(2, •U)«08 

445 


V 

9 

«.*0f ♦DS*V 

444 


If (J -N) 

21i31«21 

447 

31 

If (KN mt) 

32.33.34 

448 

32 

V 


(6,*VN •Y)/U 

444 


CO TO 3S 



470 

33 

V 


VN 

471 


CO TO 3S 



472 

34 

V 


(0$*VN •FPP(I))/(t. ♦Ff(l)) 

473 

3$ 

e 


V 

474 


D 


0$ 

475 

«! 

DS 


s(.n -sn) 

476 


U 


FRP(I) •(•R(l)*V 

477 


rppp(i) 


JV -Ul/OS 

478 


ppp(i) 


u 

474 


rp(i) 


(f(j) -Fdn/os -osacv 4u ♦u)/<*. 

480 


V 


u 

481 


j 


1 

462 


I 


I 

483 


IF <J •«) 

RlrSSRl 

484 

51 

FPPP(N) 

a 

FPPRTn.1) 

485 


FPP(^^) 

a 

B 

464 



a 

Of ADaIFPRIN*!) «B 

“87 


1 

s 

1 

488 

St 

RtTUR^^ 



484 


END 



«40 


subroutine 

ABCL (S.STEP.Xl.X2,X3,Yl,Y2.Y3.H.N,TnL) 

«4l 


CALCULATES ARC LfNGTH USING AIRST THREE OEPIVlATIvES Of X ANO X 

492 


OP 

3 

STEP 

493 


p 

9 

,5*0P 

494 


N 

9 

1 

495 


S 

e 

SO^’tYl**? ♦Yl«*2) 

496 


XX 

s 

XI ♦8TER*(X' ♦,5«STEP*X3) 

497 


YT 

s 

Yl ♦STEP*(Y2 ♦,5*STtP*Y3) 

490 


$ 

B 

S ♦SQRT(XX«*2 ♦YY**2) 

499 


XX 

B 

XI ♦P«(X2 ♦,S*P*X3) 

500 


YY 

a 

Yl *P*(Y2 ♦,5*»*Y3) 

505 


SUM 

9 

80RT(XX*»2 ♦YY**2) 

502 


SUM 

9 

8l'N*0P«2,/3, 

503 


8 

9 

SUN «8*0P/6. 

504 


DO 12 1«2< 

,M 


$05 


81 

9 

8 

506 


8 

9 

,5*(S •,5*SUM) 

507 


OP 

9 

,5*0P 

508 


P 

9 

,S*OP 

504 


XX 

9 

XI ♦PttX? ♦,5*P*X5) 

510 


YY 

9 

Y! ♦P*(Y2 ♦,5*P*Y3) 

511 


SUM 

9 

S0RT(XX**2 ♦YY««2) 

512 


N 

9 

2*N 

513 


L 

9 

N -1 

5ta 


00 14 K«1«L 


515 


P 

9 

P 40? 

514 


XX 

■ 

XI *P*(V2 ♦,5*P*X3) 

517 


YY 

9 

Yl *P*(Y2 ♦,S*P«Y3) 

r\a 

>4 

SUM 

9 

SUN «SQRT(XX«*2 «VY*«2) 

514 


SUM 

9 

SUn#dp*2,/3, 

520 


8 

9 

8 ♦SUN 

521 


ERR 

9 

8/81 •!. 

522 


IF (ABSffHR) •TOt) 2J,21il2 

523 

12 

CONTINUE 



52« 
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p 


1 


I 


( 


I 


I 


1 


c 

c 

c 


It RCTUAN 
CNO 

SUBROUTINE eSTlN (R« 10« IHAX# JNAX) 

••••••CtvCB INITIAL EITIHATE CF ROUnTIAL A 8 ZERO PERTURBATION—— 


DlHENtlON R(ID,n 
DO ttO IaipI»>AX 
00 «0 JaipJNAX 
40 RtIiJ)RO. 

RETURN 

END 

OVERLAV(JERRV,l#n 
program ONEI 
COMMON P(B1.B2) 
common XR(M),YB(«1)»CP(«I) 

COMMON thet(si)« •HETR( ei)«sT(ei)*CTrM)fMi («npR2(Et)*«3(en 

• #M«(8i),M5(en»YHP(«n.opc(6n»P(8n,M<(Bi)«8(en 
COMMON /BLOKl/ X3T 

COMMON /BL0R2/ PlpRAO 

COMMON /BLOKj/ IM4X,0UMMY(17),ML0SE 

COMMON /RL 0 M 5 / JMl» 0 Y,n 

common /BU0M6/ *0(100), YO(10C),XOP(lOO)pXORP(100)»*OPPP(JOO), 

• YOP(IOO) , YOPP(IOO), YOPPP(IOO) ,S00(100) ,I*Y,OYOXN,OYO*T 

COMMON /BLOK 7/ S*'**,SI,XM,kIm,a«,0X10*O,0X10*m, A?,* 5,X’0,*ll ,C*M, 

• DX,X10,XREF 

DIMENSION XBl (100),YBI(100),*B2(100),YM2(100) 

dimension 01 (ion,o2(ion,DJnoi),o«(ion,o5non 

IF(KtOSt,EO,0)r.O TO 10(» 

lUlMAK 

SMAXbSOO(IXY) 


Ac(3,*SM*X*OXIOXn)/2, 

B«(I,*(SMA*«A) 

0Y*1,/(Imax«1) 

xx«o, 

00 I I«1,IMAK 

8(1 )«,5*SMA*t(*x«,5)* tAY8*(X*»,S)**2) 

0XI0XaAY3,*B*(X*-,5)*«2 

F(nol,/ 0 X 10 * 

xxaxx^ox 

1 CONTINUE 

CALL INTPUl ,IHAX,S,XB,XBI ,*fi2, 1 , IXY , 30J, *0, XOP, *OPP, XOPPP ) 
call INTPLd fIMAX.S, YB,YBi, YB?, 1 , 1 X Y , 800 , YO , YOP , YOPP , YOPPP ) 
DO <1 IB1,IMAX 

AK(n*8QRT(X8?(n»*2*YB2(n**2) 

IF'(XBl(n,LE,l,)G0 TO 2 
NRITf ( 6 ,RR 0 l)>ai(n,I 
XBHDmI, 

2 CONTINUE 

if(abs(ybi(1)),le.i,)go to J 

NR1TE(6»RR02)YB1 (I),I 
YBl(l)«8:GN(t,,Y8l (I)) 

S CONTINUE 

THeTX«Sir,N(ACOS(*01(I)),YBl (in 
TMCTY«*SIN(YBl(I)) 

THFT(I)*,5«(TMfTX4TMETY) 

4 continue 
THET(n«. 5 «Pl 
THETnMA*)».,'i«PI 
RETURN Ak 


52^ 

52B 

52T 

S2B 

529 

530 

531 

532 

533 

534 

535 

536 

537 
536 

539 

540 

541 

542 

543 

544 
54^ 
54fc 
547 
546 

549 

550 

551 

552 

553 

554 

555 

556 

557 
556 
a59 

560 

561 

562 

563 

564 

565 

566 

567 
566 

569 

570 

571 

572 

573 

574 

575 

576 

577 

578 

579 

560 

561 

562 

563 

564 
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100 CONTINUE S8S 

CAU MOHTANaxV,XO,YO,800,XOP,XOPP,XOPPP,VOP.VOPP,VOPPP,SHT,>tMT 586 

• »yHT,108HT) 587 

XSTaXOnxY) 588 

aiaSHT 589 

XlOaXMT 590 

CALL 8ETuP0(lMAX«Sl«XH,X|M,A4*DXlDX0,DXl0XM,A2,A5»CXM,0X.Xt0) 591 

OXal./HMAX-n 592 

XXaO. 593 

00 101 Jal.I^AX 599 

call SDRIV0(XX»SS«DXI0X«XIM«49,DXI0X0«DXI0XM, A2, A3«CXN) 595 

SdJaSS 5«6 

F(I)*1 ,/nxiOX 597 

IF(S(I).L£.Snilol 598 

XXsXX+DX 599 

101 CONTINUE fcOO 

CALL lNrPL(l.Il»S.XB,X81,Xe2,l,lXY»SOO,XO,XOP,XOPP»XOpPP) 601 

call lNTPLn.l1»S,VB,YPl.V82,l»lXY,SOO,YO,VOP»YOPP.VOPPP) 602 

00 10« 1=1,11 603 

AKn'8SQRT(xe2(n=*2 + Ye,*:(I)**2) 60tt 

IF(X01(1),LE,1,)OO Tn l<2 60S 

WRITE(6,9901)XBl(n,l P06 

XBKDsl, 607 

102 CONTINUE 608 

lF(A6S(YBUin.LE,l,)G0 TO 103 609 

WRlT£(6,9902)Y81(i)»l 610 

YBl(I)«SIGN(l,,YBim) 611 

103 continue 612 

THETX«S16N(AC0S(XBl(n),YBl(n) 613 

ThETYaASlN(YBim) 619 

THCTn)».5*(THETX*Tht:TY) 615 

10« CONTINUE 616 

TH£Tm = ,5*PI 617 

ii (I09HT,LT,IXY)G0 TO 105 616 

12=11 619 

60 TO 111 620 

105 CONTINUE 621 

10T=108HTtl 622 

01(10T)*XHT 623 

02n0T)=VHT 629 

00 107 I=IOT,IXY 625 

Dl(l*l)8X0fl) 626 

02(lAl)=Y0n) 627 

107 continue 628 

IXYP1=IXY+1 629 

CALL 8PLlFnOT,lXYPl,01,D2,O3,Dfl,O5,l»0,,l,OYOX1 ,IN0) 630 

1KAXK1=1NAX»1 631 

UP1«IH1 632 

00 108 IeIlPl,IHAXBl 633 

D3=S(1)«SHT 639 

XOmaXKTAOS 635 

IF(XB(I).GT,XO(1XY))GO TO 108 636 

I2«I 637 

108 CONTINUE 636 

CALL lNTPL(IlPl,l2#XB,Y8,yBl,Ye2,IOT,lXYPl,01 ,02,03,09,05) 639 

on no i=iipni2 69o 

AK(l)a«Y82n) 691 

THETB(1)sATAN(YB1(1)) 692 

7METtn=0, 693 

110 continue 699 
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111 CONTINUE 
12P1«1*41 

DO 112 IaI2Pl»lNAX 
08a8(I)»SHT 
XB(I)«KHT408 
YB(T)«tO(lXY) 

AK(I)aO| 

TMET<mO. 

TMETBdJaO, 

112 continue 

RETURN 

R901 FDRMAT(* X61»*Elt>,«* AT Ia*13) 
9902 FORMAT!* YPl (n«*EU,e* AT I»*13) 


END 

SUBROUTiNf InTPL(mI»nI,SI*FI#P1P»F1PP»H»N,S»F»FP»FPP»FPPP) 
interpolation U81NG TAYLOR SERIES • JAMESOn 

DIMENSION SHn»FI(n*FlP(l)#FIPP(l)»8(l)»F(l)»FPn)#FPP(n 
»FPPP(1) 



K 

c lAflStN »M) 


K 

■ (N »M)/K 


I 

B M 


MIN 

a MI 


NIN 

s NI 


0 

B S(N) *8(6} 


IF (0*! 

8i(Ni) •si(Min) n 

11 

MIN 

a NI 


NIN 

B Ml 

13 

K1 

B IAPS(N1N -MIN) 


IF (KI) 

?1,21,1S 

15 

KI 

B (NIN «MIN)/K1 

21 

11 

a MIN aKl 

31 

II 

a 11 *K1 


88 

a Sldl) 

33 

I 

a I *K 


IF d 

-N) 35,37,35 

55 

IF (D*(3d) -ssn 33,33,37 

37 

continue 


1 

a I »K 


ss 

B 88 -8(1) 


11»13>13 


FIPPd DeFPPd )*SS*FPPP(I) 
FiPdnsFP(n+s8*(FPPd)*ss*FPPP(n*,5) 

FI dnaF(n*SS*(FPd)*,5*SS*(FPPd )*SS*FPPPd )/3, ) ) 

IF dl -NIN) 51»U1»51 
41 RETLRN 
END 

SUBROUTINE SETUPO(IMAX,S1»XM,XIm# A 4, A1 ,BB» A2, A3*CXM,0X, XlO) 

XIM8SJ»XM*X10 

OXal, /(lMAX-1) 

Clax:M/CXM*Al 

C2aBB«Al 

C3b2,*CXm*PB/(1 .-CXM) 

X2aCXM«*2 

X4bX2«*2 

X6ax<i*x2 

A2a(70,*Cl«22,*C2*2t*CS)/l6,/X2 

A3a(»eu,*Cl*36,*C2»Rt*C3)/16,/'X4 

ARtCSO.aCl-lR.aCE+a.tCSl/lOi/XB 

RETURN 

END 

SUBROUIINE 80R1VO ( XX , S» OX IDX , X I H , A4 , A I , (3B, A2 , A3» CXM) 


645 

646 

647 
64# 

649 

650 

651 

652 

653 

654 

655 

656 

657 
65B 

659 
6^0 
66 1 
662 

663 

664 

665 

666 
o67 
668 

669 

670 

671 

672 

673 

674 

675 

676 

677 
676 
679 

660 
66 ] 
662 

663 

664 

665 

666 

667 

668 
669 

690 

691 

692 

693 

694 

695 

696 

697 

698 

699 

700 

701 

702 

703 

704 
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»CAUCUL*T?$ 8 and OXIDX AS FUNCTIONS OF X 


IF (XX.CT.CXM) 60 TO 10 

X2aXX**2 

X4«X2*«2 

X6>X««X2 

TlaA2*X2 

T2bA3*X0 

T3aA4*XA 

8axX»(AUTltT2tT3) 

DXl0XeAl*34aTl+5,*T2*7,*'’3 

RETURN 

10 IF(*HS(XX-l , ) ,LE, 1 ,t-06) GO TO 20 
TSXX-CXH 

TJai,/n,-T/(l,*CXM)) 
SsxlM^BB*T*Tl 
DXIDXaBB*Tla*2 
RETURN 
20 SS 1 .E 430 

DXI0X81.E+30 
RETURN 
END 

SUBROUTINE 


n t 1 

ENO 

SUBROUTINE KORTAN(N,X»Y#S»XP»XPP»XPFR»YP»TPP»TPPP|SHT,XHT,VHt»lH) 
dSSJSsIJn xm;Y(n:s(n.xP(;),xPP(n,xpPP(i),YP(n.vpp(i).YPPP(i) 
DIMENSION 8M(n I XM( 1 ) » YH( 1) *0UM(l ) 

DO I I»1»N 

IF(YP(I).LT,0,)GO TO 2 
IHal 

1 CONTINUE 

• continue 

8HTaS(IH)tYP(I^^)/(YPnH)*TP(I^+l))*(S(IHil)«S(lH)) 

SH(l)aSHT 

CALL INTPLCI * 1 ,SH,XW»OUM,DU»' , 1 , N,. S . X , XP, XPP, XPPP) 

CALL INTPL(1»1 ,SH,YH,DUM»0UM,1,n,S,Y» YP, Y pp^Yppp) 

XHTBXH(l) 

VHTaVHt 1 1 

RETURN 

END 

LVERLAY(JERRY, 1 ,2) 

PROGRAM ONE? 

COMMON P(6J#82) 

COMMON xp(en»YB(8n»cp(6i) 

COMMON THET (B1 ),THET('(81)*ST(ei)>CT(6P»Ml(6l)»*2(en»'*3(61) 

* ,w<i{en.ws(6i),YBP(eu»opo(8n»F(8j).AK(6n,s(8n.LS(nn,fM{6n 
COMMON AN(8n,G(Rn»GM(81),CB(6n#0(8n»Xl(81.1.X2{en»M(8l),HR(81) 

l,HRP(ei) .HRH(61 ) ,HRMM(6n 
COMMON XS (LOO) ,Y5 (POO) 

COMMON ID, ANMAX,DNDY0» YMAX,CD,NHSQ, JSMP 

common /BLOK2/PI»RAD/BlOK3/ImaXiJmaX,C2»RE1»DPM,IOP*J 

IDP.RPM, JR, JR,NS,GM102» AOSQ,DXSOf OXnY,OYSO»PX?,OY2,KLOSE 

COMMON /BLOKS/ JM) ,OY, 1 1 , JSUP, JSON* OF 3 , 1 SUB , 15UP»SUMRP 

common /BLOKT/ 6MAX,SliXM,XIM,Aa,DXlUXO»OXlUXM,A2,A3»XlO,XIl»CXM, 

* OX 

COMMON/0LON6/ALF 

CALL NTRAMF (AN, ANMAX, JHAX,ONDYO»OY,G»GH, ALF) 

CALL k1m?(TMFT,THETB, YB,YBP,"1 ,K2,w3, Kii,KS,ST,CT»C( JHAX) ,ONOYO 
»DY» n » 1MAX,KL03E»ALF) 


RETURN 

end 


705 

708 
707 
70e 

709 

710 

711 

712 

713 
7l« 

715 

716 

717 

718 

719 

720 

721 

722 

723 
7?R 

725 

726 

727 
7?e 

729 

730 
751 
732 
755 
750 

755 

756 

737 

738 

739 

700 

701 

702 
75 3 
7o« 

705 

706 

707 

708 

709 

750 

751 

752 

753 
750 

755 

756 

757 

758 

759 

760 

761 

762 

763 
760 
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SUBROUTINE NTRANF (X,XMAX, JHiX,ONOVO»OV.O|6H.Atr) 76$ 

766 

•••«-«COMFUTES STRETCHING OF NORMAL COORDINATE-— TaT 

768 

dimension X(t)» G(l)» 0H(1) 764 

B=l,/0N0Y0 770 

IF tXMAX,GE,l,C*06) CO TO 10 771 

AaB»l,/XMAX 772 

K«0 771 

CO TO 20 770 

10 K*1 775 

AsB 776 

20 Oyal,/(jMAX-l) 777 

00 50 JeliJMAX 778 

ZETABl,-(J-n*OY 774 

IF (JBK.FO.n GO TO 10 780 

AA«(1,-ZETA)**ALF 781 

X(J)a2ETA/e/AA 782 

GO TO 40 781 

10 GdlaO, 764 

GO TO SO 785 

40 C(J)=8*AA4a,-?ETA)/(l,-(l.-ALF)*ZETA) 786 

IF (J.EO.n GO TO 50 787 

CH(J-l) = ,5*(C(J)*G(J-n) 786 

50 CONTINUE 784 

AA«(1,+DV)**ALF 740 

GJPlaB*AA*(l,*OY)/(l .♦(1,-ALF)*DV) 741 

GH(JHAX)8,5*(GJPUG(JmaX)) 742 

RETURN 745 

END 74u 

subroutine MH2(THET,THETfl,YH, YBP,*.l.*<2,v.3,lNU,t»5»ST,CT,G,0^DYC 74$ 

* .CY,Il*IHAXiKLOSe,ALF) 746 

747 

-•calculates Y8p(n,*.i(n,M2(n,sT(n»cT(n 748 

744 

OlNENSION TBETtn, TUETBU), YR(n. YBP(U. ( 1) * ^2(1 )» "KD. n 6 no 
14(1), M5(l), STH), CT(1) 641 

DO 10 lai,lMAX 802 

8T(naSlN(THET(l)) 603 

10 CT(n=COS(TMETfl)) 804 

00 20 l«l,tl 60S 

YBPCUsO, 806 

Wl(I)al, 807 

W2(l)a0, 608 

W3(I)sl, 604 

W4(1)S0, 610 

W5n)s2,aOY/C 611 

20 THETB(l)»THETm 612 

IF (KU03E,EQ,n RETURN 613 

IlPlall+l 814 

00 30 1411F1,U'AX 615 

YDP(l)*TAN(TH£7e(n) 616 

Ml(I)«l.*TBP(n**2 617 

w2(nBYBP(n/wi(n eia 

Ylo-OY 614 

1F(YO(I),CE.ONOYO)GO to 25 620 

YBOAsYDO )/DSDYO 621 

YA»-YB0A-ALF*YB0A'»k?«iALF*( 3,«*LF-1 , )/a,*YB0A<"*3 C?2 

* •ALF*(l6,*Ai.F**2-i2,*ALF*2,)/6,*'Yb0A**4 623 

IF(.'HS(YA),GT,OY)GO TO 25 624 
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VI«0Y4?,*Y* 

25 COMTINUC 
YlOOY^Vl/OY 

NUl)aV} 00 v **2 

‘<S(nB«Yt*(l,*Y100Y)/C 
30 CONTINUE 
RETURN 
ENO 

0VER'.AY(JERRY,1,3) 

PROGkAM ones 

c 

c I'lON OP POTENTIAL eon, PY COLUMN REU A X A T I ON- * •••• 

c 

CONNON • (6l,fl2) 

Cr-NMOM x,'(fli),vB(en.CP(6l) 

cohmun THET(pn,THETfi(en,ST(fln,cT(0n,Niten»"?(fln,''S(Hn 
* »'‘«(8i)r»'S(en» tbp'bi) » opo(8i) ,F («n » AKf PI ) ,s(Pi ).LS(Pn.p''(en 

l,HRP(fln,HRH(8n,HRMM(Bn J / V i; 

COMMON /BLOf'S/ I'^AX,J^'AX,C^,RFl,0^'H,I0p,J0P,Rp^.,IH,JR,^S,GMlO^ 

J, A0Sri,0XS0,0XDy,nvSQ.DX2,0Y2,KL0SE 

COMMON /BLOKE/ JMl,OY*n, JSUP#J80N,0F3,18 uB,ISUP,SLimrp 

SUMRPao, 

0F1»1./RF1 

J0a2AJ3UP 

IF (JO.EQ.Z) JC32+JS0N 

JPlcJHAX^i 

JsJMAX 

DO lO 

10 0P0(I)s3T(n 

IF (Kucse.Eo.n CO to so 
I lPlall+l 

IMAXMlsl-AX.l 
DO 2o I»I IPl » I«AXMi 

20 0P0( J)3W2(i)*( t,*A ( ’)<*0X?«(P(U1,J)-Pn.i,j))) 

DPIs3,*P(IMAX,JMaX)»U,*P(IMAX*1,JMAX)*P(IMAX-2,JMAX) 
0P0(IMAX)sw2(IM4x)«( I ,+F( lMAX)*0X2*0Pl ) 

30 CONTINUE 
C 

c .start a cycle at imCSTAC, PT, 

c 

OPM30 , 

RPMaOPM 

I=l 

jsjo 

NSoQ 

KS»0 

IIMls*i.i 

A6«0, 

Blso, 

BaeO, 

GO TO 230 
«0 IMMbI»2 

IF (I, EG. 2) IMMt2 
C 

C -r,.., .COMPUTE QUANTITIES DEPENDING ON I ALONE 

C 

FDefnT*0X2 

FOl5F(4)*0XSQ 


625 
Q26 
627 

626 
626 

630 

631 
652 

633 

634 

635 

636 

637 
e:" 

839 

640 

641 

842 

843 
6aa 
bas 

8C6 

day 

eae 

ea9 

650 

asi 

652 

655 

65a 

655 

656 

657 
656 

659 

660 
661 
662 
665 
66U 

665 

666 
667 
666 

669 

670 

671 
872 
673 
67a 

675 

676 

677 

678 

679 

680 
661 
662 
663 
66a 


30 


F02«f m*OXOV 

KS«0 

cc«o. 

DOxCC 


-•-•••COmFUTC QUANTlTieS OEPENOlNf, ON 1 *NO J — 

50 CONTINUE 

HRm(J)sMR(J) 

. MR(J)8MRP(J) 

IF n.GT.HNl) on TO 60 
HRP(J)sJ,/(l,*AK(I*n*AN(J)) 

GO TO 70 
60 HRP(J)*1, 

ASaO, 

70 CONTINUE 

S»aP01*HR( J) 
saaF02*C(J)*HR(J) 

S3aG(J)*0VS0 

S4aG(J)*DVi! 

S5^F0*HR(J) 

TInsF(I-i)*mRH(J> 

FHMa,S*(HF*Tl^') 

AKHsAK(n*H»(J) 

RRtl,/(YB(n*AN(J)*CT(l)) 

• --..•CONPUTf PHI«DIFf ERENCES FO® VEUOCITY COMPONENTS, NOTE 

increasing j means decreasing norhalczeta or n) coordinate, »-•- 


0P1SP(IY1»J)-X1(J) 

opjsp(i.j-j)-p(i»J+n , _ 

0Pjj8GH(j»i)*(p(i*J“n*'p(i»j5)»GH(j)*(P(i*j)»P(*<j+n) 
UsCTd )f0PI*S5-YHP(n*0PJ*SU 
Vs-ST(l)+DPJ*S« 

VBaVYBPCl )*U 
LSI 

IF (V,ET,0,) La-t 
IF (J,6Q,JNAX) Ua.l 
TsU 

UUaUAU 


VVbVaV 

OQsUU«VV 

AAaA0SQ«6H10E«0Q 
ARal ,/AA 
T4B1 ,«UU*AR 
UVsu*V 
UVARcuv*AR 

A<ia(AKH*TU + nR*CT(n)*8« 

IF (I.GT.IU on TO 80 
A5a (a,*AKH»UVAH»RR*ST (!) )*SS 
80 CONTINUE 

rHe,5*(HFPF(I*n<*H«P(J)) 

Dpijsp(i4i » J«n-P(i*i»J+n4p(i-iiJ^i)«pn-i » J»n 

opiieFH*(P(i4i , jn"FHM«(P(i, j)-p(i«i»j)) 


BIBO, 

OflRO, 

A6m0, 

ir (0Q,1.T,.''A) GO TO 120 




665 

666 
667 
666 
669 

690 

691 

692 

693 

694 
895 

696 

697 
696 
699 

900 

901 
9C2 

903 

904 

905 

906 

907 
906 
969 

910 

911 

912 
9)3 

914 

915 
9)6 
9)7 

916 

919 

920 

921 

922 

923 

924 

925 

926 

927 
926 

929 

930 

931 

932 

933 

934 

935 

936 

937 
936 

939 

940 

941 

942 

943 

944 
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aa.MMaRAric u 4pn ftPcniun httfa fhr ftupPDQniuTr n 

945 



947 


KSaKS^t 

948 


0Rb1,/QQ 

949 


8US1*0H 

950 


8a«S2*QR 

951 


$3sS3*OR 

952 


AUUbUU«S1 

953 


AVVaVV*Sl 

954 


uBMU«rRP(i)*v 

955 


BUVNbS2*V*UB 

V5t) 


cuusub*ub*$3 

957 


CVV*VB*VH*S5 

956 


IF (QF3.LE,1 ,E«06) GO TO 91 

959 


FAC«SQRT ( ARS( 1 .■00/AA) ) 

960 


BIa«QF3*FAC*ABS(VR)AG(J)*DXDY 

961 


8«b0F3*F4C*U*F0I 

962 


AbsB<l*(P(l*lfJ)*Xl(J)) 

963 


continue 

964 


OPNNaAVVtOPI 1«BUVN*DPIJ+CUU*DPJJ 

965 


KMs(J*J«ifL)/2 

966 



967 


IF (JM.GT.n GO TO 100 

966 


JRM»1 

969 


KMMsl 

970 


GO TO no 

‘>71 

too 

KMMSKM+L 

972 



973 

no 

CONTINUE 

97(1 


FHHMa,5»(TIM»rnNM)*HRMM(J)) 

975 


0PIl6FHM*(P(I,J)-P(I-l,jn-FHM-4(Pn,l,J5«X2(J)) 

976 


DPIJaPd, (JN) 

977 


BUVSb«<i,4S2*T-»U*VB 

978 


A2SaGH(KMN)+GH(KM) 

979 


0PjJtGH(K«H)*P{:,JMM)«A2S*P(I. JH)4GH(KW)*P(1, J) 

980 


DPSS»AUU*CPI I'flUWSaOPIJ+CVV^OPJJ 

981 


AlSsi ,«QQ*&R 

982 


B2s,5«A1S*(A2S*CVV»BUVS) 

9fli 


B3aRlffl2 

984 


A»( l,»T)*fa3«CUU*GH(J-J)"Aa 

985 


Ca(l,4T)*P3«CUl'*GH(J)*Aa 

986 


Ba-A»CfAlS*(HuvS»2,*AUU*FHH)*AVV*FHM4BU 

987 


RP8AlS*0P8S*l)PNN4Aa«0PJ+A5*DPl**b 

988 


ARP84BS(HP>Ab) 

989 


8UMBPaSU'’‘RPY*«P 

990 


IF (ARP. lE, PPM) CO TO 140 

991 


ISUPaj 

992 


ISUBcO 

993 


IRal 

994 


JRaJ 

995 


RPM8ARP 

996 


.GO TO 140 

997 

120 

AlsTc*S 1 

998 


A2ano*vnp(n»uvARi*s2 

999 


A3n(*> 1 (I ' -AR*VB**2)*S3 

ior»o 


HPaAl*OPn-A2*OPIJ + A3*OPJJ»A4AOPj4A5*OPl 

)001 


ARPeABS(RPI 

1002 


8UNRP»SU”RPaARP 

1005 


IF (ARP.lE.RPm) GO TP 130 

1004 
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isueal toos 

I3U»«0 t006 

IR«t 1007 

JR«J 1006 

RPMIARF 1000 

1010 

••.•••COHPUTC TRIDlAGONAt COEPFS— 1011 

1012 

130 A«*A3*GH(J«1).A0 1013 

Ca«A3*GH(J)«AA 1014 

Ba*A«C*QFl«A1*(FH4FHM) 1015 

140 CONTINUE 1016 

e»l,/(B«A*CC) 1017 

CC«B*C 1016 

D0«B*(NP«A*D0) 1010 

IF (J.EQ.JNAX) CD TO ISO 1020 

CB(J)«CC 10/« 

0(J)«0D 1022 

J«J«1 1023 

60 TO 50 1024 

150 0P*OO 1025 

IF (ABStOP).LE,OPN) GO TO 160 1026 

CPH«A0S(OP) 1027 

lOPal 1026 

JOPaJ 1024 

160 X2U)sXl(J) 1030 

XUJ)«P(I«J) 1031 

PU»J1«P(I.J)40P 1032 

00 190 JJaJO,J“l 1053 

JaJ«i 1034 

OPaO(J)*CP(J)*OP 1055 

IF (ABS(0P),LE,0PN) 60 TO 180 1036 

OPMaAHS(OP) 1037 

lOPal 1036 

JOPaJ 1039 

180 X2(J)aXl(J) 1040 

Xl(J)aP(t,J) 1041 

190 P(l»J)aP(l,J)*pp 1042 

Jajo 1043 

LSmaKS 1044 

NBaNSAKS 104S 

p(i»jMAX4i)8W3m*P(i.jM4x*i)4N4(n*pn»JF4x)*>«5m*DP0(n io46 

1047 

••••••CHECK I FOR ENO OF CYCtE.IF BOtY IS CLCSfcO, lalMAX IS SYhN£TR 1048 

AXIS, IF BODY IS OPEN, lalMAX IS EITHER NOT COHPuTf 0 (SUBSONIC FREE 1049 

STREAM, P(IHAX,. 1)50,1 OR EXTRAPOLATEO(SIJPERSONIC FREE STREAM)—— io50 

1051 

IF (1,E0,IMAX) RETURN 1052 

lal*l 1053 

IF (I,EQ,IMAX) GO TO 200 1054 

GO TO 40 1055 

200 IF (KlOSF.EO.l) CO TO 230 1056 

IF (JS0N*J8U?,LT.1) 60 TO 220 1057 

00 210 JaJO,JPl 1058 

210 P(I,J)a3,ACP(I-l,J)«P(l»2,J))*P01»J) 1059 

220 RETURN 1060 

1061 

.••—SPECIAL EONS FOR SYMMETRY AXIS»I = 1 OR Imax— — 1062 

1063 

230 CCaO, 1064 
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DD«CC 

$l«2,*DXS0«Fn)**2 
S3«AK(n*0Y2 
IF (l.CO.n GO TO 250 
00 240 JJa2,JHAX 
240 HR(JJ)«h«F(JJ) 

L*»l 

TAaO. 

TCaj. 

. IF (I.EQ.IMAX) GO TO 270 
250 1M82 
I «» 

TAbi, 

TC«0, 

DO 2«0 Jjs2,jKix 
HR(JJ)al ,/n ,♦4K(n•A^(JJ} ) 

HRP(JJ)sl ,/( J .♦ aK (2)**N(JJ) ) 

HRH(JJ)*MRP(JJ) 

260 Xl(JJ)8p(2,jJ) 

270 OPllsPdM# 

DP.i»p(i,j*i)»pn» J+n 
V»*STn)+0V2*C(J)*DPJ 
VV8V«V 

AAbA0S0-GH10?*VV 

COMPUTE COEFFS OF OIFF £0, AT SYMMETRY AXIS------ 

AJr2, 

A3=83*41*G(J) 

At*Ai«Sl*HR(J) 

A2a(l,-VV/4A)*C(J)*0YSQ 

BlaO, 

KSlJP*0 

KSUBbO 

IF { J , EQ , J‘’l . AMO . I . tO . lMAX ) GO TO 290 
IF (VV.GE.AA) go to 300 

290 OP;jeGH(J-l )*P(I, J-n-(GH(Jr.j )4GH(J))*P(I, J)fGH(J)«P(l, j*n 
GO TO 310 

300 continue 

IF (QF3,LE,ltE-'’6) 60 TO 301 
FAC«80RTCABSn-VV/AAn 
0l«**BS(V)*F*C*G( J)*2,*OXOr*OF3 

301 CONTINUE 
KSUPal 
KSRKS^l 

KMb{J*j-1.U)/2 

KMNaKM-l, 

JM8J.U 

jMHaUHaU 

A2S«r.H(KMH)tGH(KH) 

0PJJ«GH(KMM)*P(I,JKM).A2S*Pn , JM)4GH(KM)*P( J, J) 

BaA2S*A?«Bl 

AaTA*B*A3 

CaTC*B*A3 

1 

GO TO 320 

310 As.a?*GHC J» i )«&3 
C»«A?#GHU)^A3 
Bo*A*CtOn*Al 


106 % 

1066 

1067 

1066 

1069 

lOTO 

1071 

1072 

1073 
1076 

1075 

1076 

1077 
1076 
1079 

1 o«o 
toftj 
10*2 
1063 
1066 
10*5 
10*6 
10*7 
1 0*6 
1 0 r 9 
10^0 
1 09 J 
109 ? 

1 093 
1096 

1095 

1096 

1097 
1096 

1099 

noo 

! mi 
1102 
1103 

1100 

1105 

1106 
1107 

t ma 

1109 

mo 

mi 

m2 

ni3 

nu 

ms 

1 U6 
1117 
m A 

1119 

1120 
1121 
1122 
1123 
1129 
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K$UBb1 

320 RP«A1«OPII«A2*OPJJ4A3«OPJ 
ARPBA88(PP) 

SU'*RPbSUmRP«ARP 

IF (ARP,LT,RPH» 00 TO 330 

RPMaABS(RP) 

ISUBaKSue 

ISUPBKSUP 

IR«l 

JRBJ 

330 CONTINUE 

B«1./(B*A*CC) 

CCbB*C 

O0a6*(RP«A*O0) 

IF (J.TO.JMAX) 60 TO ISO 

CB(J)«CC 

D(J)B0n 

J*J41 

GO TO 270 

END 

OVlRLAV(JERRY«l»a) 

program ONEO 

COMMON P(81«B2) 

COMMON xB(8n#Y3(Bn»cr(en 

COMMON THETIBn.THETBfBnfSTIBn.CKBn^wKBn.MPfBnfNKBt) 

* ,M4(8n,M5(8n#Yep(BU»opo(8n,F(Bn»AK(8i).s(ei),LS(en,FM(en 
COMMON AN(6l)»G{8l)»GH(ei)#CB(8U»0{8l)*XI(8l)»X2(8n«N(8t),MR(8n 

I,MRP(81)*MRM(8n»HRMM(8l) 

COMMON XS(«00) ,VS(000) 

COMMON 10, ANM*X,ONOYfl,YMAX,CO»RMSO,JSM?,TSP 
COMMON /BLOKl/TMiX, JM*X,C2,RF1 ,0PM; lOP. J 
;OP,RPM, IR, JR,NS,r,Mi02» AOSO,OXSG,nxOY,OYSO,OX2,OY?,KLOSE 
COMMON /BLOKfl/ GMSfl, 60GM1 , T0CM80, CPO, MST AR 

COMMON /ML0K7/ SMAX,SI,XM,XIM, A«,OX10XO,0*IOXM,A2, ASjXlOjXll.tXI, 

* OX,XtO,XREF 

call CPBOOY (P,F,Ml,YBP,0PO,CT,LS,CP,fM,l0,lMAX,JMAX,CMlO2 
« • AOSO,OX2,KLOSE,GMSO,GOGM1 ,T06NS0) 

KRIT£(6,570) 

WRITE(6, 580) (I, 3 (I),XBm,YBm,CP{n,FM(n,I*l,IMAX) 

CALL DRAG(r.P,Yfl,TMET,THET6,r ,RMSQ,1MAX,0X) 
call CPPLOT (S,X«, yb,CP,IMAX,CPO,KSTAM) 

RETURN 

570 format nHl,8X,lMI,6X,2HSB,8X,2MXB,8X,2KYB,ex,2HCP,8X,lHM/) 

580 format (no,3F10,3,2F10.S) 

END 

SUBROUTINE CPBOOY (P,F,WJ , YBP,OPO,CT,LS,CP,Fm, 10, IMAX, JM4X,Cmi02 

* ,a0$0,0X2,KLO$E,GMSG,COGM1,T0GMSQ) 

c 

c —•••COMPUTES SURFACE PRESSURE COEFFICIENT ANO MACH NO,-*— 

C 

DIMENSION P(10,l), F(U, Wl(l), YBP(J), OPOO), CT(1), LSd), CP(l 
1), FM(l) 

C 

JBJMAX 

OaO, 

00 60 I«i,IMAX 

IP-lAl 

IPPBI42 

IF (I,EO,lMAX-n IPPel 
IM«1*1 


125 

126 
127 
12B 
12S» 

130 

131 
<32 
133 
13« 

135 

136 

137 

138 

139 
IRO 
1«1 
H2 
103 

14UI 

105 

106 

107 

108 
100 

150 

151 

152 

153 
150 

155 

156 

157 

158 
156 
160 
161 
162 
163 
160 

165 

166 

167 

168 

169 

170 

171 

172 

173 
170 

175 

176 

177 

178 

179 
IPO 
181 
182 
185 

180 
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lies 

irCX.CC.DGO TO «0 1186 

IF (I.EO.Imax) CO TO so 1187 

IF n,tO,2) 1MM«2 1188 

0J8F{1P,J)-P(IM,J) 1106 

GO TO SO 1190 

20 continue 1191 

0J95,«P(I,J)-«.*P(IM,J)4P(Imm,j) 1192 

C .compute surface VELOCITV.-.— H9J 

30 u3CTm40j*F(n*D*?.TBP(r)*t)P0(n n9« 

Q»S8BT(**1(I))*U 1195 

40 00>0*0 1196 

AAaAOS0«GM102*O0 1197 

C 1198 

c -•-•••SURFACE mACn no, •-•-•- 1199 

FM(1 )aSQRT (08/AA) 1209 

C 1201 

C -..—pressure RATIO——. 1202 

P0PINF«( 1,^CMS0*( 1,-00) )**C0GM1 1201 

C 1204 

C -. ——PRESSURE COEFF, — — 1205 

CP(I)«T0GMS0*(PQPInF-1 ,) 1206 

GO TO 60 1207 

C 1 208 

C —..IF lelMAX IS NOT A SYMMETRY AXIS, USE BCKwo OIFF F0» OPSb---- 1209 

SO IF (KtOSE.EQ.O) CO TO 20 1210 

C 1211 

C — ....lalHAX IS A SYi'METRY AXIS------ 1212 

0»0, 121? 

GO TO 40 1214 

60 CONTINUE 1215 

return 1216 

End 1217 

SUBROUTINE ORAr.(CP,R,THET,THFTH,F,RMSO,lHAX,DX) |2ie 

C 1219 

C --- — COMPUTES tlRAG coefficient BY INTEGRATION OF SURFACE PWESSu»fc 1220 

C 1221 

dimension CP{l),Rm»THETm,THETBn)»F(l) 1222 

C 1221 

C 1224 

c — — trapezoidal integration—— 1225 

C 1226 

8UM«0, 1227 

DO 10 I«2,IMAX 1228 

10 8uMosuM.(CP(n*R(n+cP(i-n*Rn-U)*(R(i)»R(i-n) 1229 

COTRAPaSUM/RMSQ 1230 

MRITE(6,540)CDTRAP 1211 

IF (MOO(IMAX,2),NE,0) GO TO 15 1232 

WRlTE(6,990n 1253 

RETURN 1234 

15 continue l£3S 

C 1236 

c -.---SIMPSON integration (ONLY IF l^AX ODD) ------ 1237 

C 1238 

SUMbO, 1239 

IMAXMlaMAX.l 1240 

DO 20 lo2,IMAXMl,2 1241 

20 SUMa8UM»CP(I )»fl( I )*SIN(TMETB( I) )/{»^ ' I) *COS ( TMt T ( I) -THE T0 ( J ) ) ) 1242 

8UMo2,*SUM 1243 

1MAXm2pI«ax-2 1244 
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\ 

I 


I 


00 ]0 I«S>1M*XH2,2 

JO 8UM«8UH*CP(!)*«(!)*81N(TMETB(n)/(» (l)*C08(THET(I)«THETBn))) 
C08lMPa«,tOX*SUM/(J,*RMSO) 

MRlTe(8«5SO)C08IHP 

RCTUBN 

540 F0RH4T (////47H OR*C COEFFICIENT OY TR4PE2010AI. INTEGRA T ION* ,F8, 
15) 

550 FORMAT {//4JH DRAG COEFFICIENT BY SIMPSON INTEGRATION*, FR, 5) 

5401 FOKMAT(/* NO DRAG BY 81HP80N INTEGRATION BECAUSE IHA« IS EVEN*) 

END 

SUBROUTINE CRPLOT (S»XB« YB«CP> Imax,CPO>KSTaR) 

DIMENSION S(l)» XB(1)» YB(n, CP(l), K00£(4)» LINEOOO) 

DATA KOOE/lH , JMA,)M0,1MA/ 

NRITC(6.S0) 

DO 10 L*l»100 
10 HNE(L)*KOOE(n 
LINC(MSTAM)«M0DE(4) 

DO 40 I*1»IMAX 
Xa4,5AJO,*(CPO-CP(n) 

IF(K,GT, 100) GO TO 30 
LlNE(K)fm00E(2) 

30 wRITEt6»70) l,i;B(I),Yem,CP(I)»LINE 
UNE(K)*kliOE(1) 

40 tlNE(KSTAR)*KO0E(«) 

RETURN 

C 

C 

50 format (3JM1PUOT OF CR AT UNEQUAL t NCREnf NT S///3X , 1 wl , 5X , ?H*B 

♦ »SX»2HYR»6X.?HCP//) 

70 format n«,2F10,3,F8,4» lOOAl) 

END 

0VERUAY(JERPY,1,5) 

PROGRAM ONES 
common P{61 ,62) 

common x0(8n,YB(ai),cP(8n 

common THtT(Rn,THFTB(en,STC8l),CT{6l),»)Cei).M?(Pl),M3(fn 
A ,MU(6n,MS(6n,YPP(en,0P0(6n,F(Bn,*K{en,s(6i),iS(6i),FM(Bn 
common AN(Bn,r,(Bn,GH(6i),ci(6i),D(6i),xt(6n,x2(fii),M(Ai),HR(en 

1»HRP(6J) ,URM(ai),HRMM(ej) 

common XS(«00) ,YS(tOO) 

COMMON 10,ANMAX,nNI)Yn,YMAX,CO,RMSO,J8KP 
common /BLOt(2/PI ,RAD 

COMMON /BLOM3/ Imax, JmAX,C2,RF1,DPm, IDP, JOP. flPM, IR, JR,N8 ,GmiQ 2 
1, AOSR,OXSO,OXDY,OYSij,OX2,DY2,KLOSE 
COMMON /BLOKS/ JMl,DY,Il,JSUP, J80N,QF3,ISUB,ISUP,SUMRP 
COMMON /RLOK7/ SMAX,Sl,XM,XlM,A4,0XI0X0,DXl0XM,A2,A3,Xl0,Xli,CXM, 

• OX 

common /blokr/ n 

CALL mcmarT (P,AK,an,F,G,YBP,0FO,ST,CT,LS,m,ID,JSKP 

* fIMAX, J mAX,GH102, AOSQ,DX2,OY2,klOSE, 1 1 ) 

CALL SONlIn (P,f ,8 T,CT,XB,YB,AK,Fm,yBP,0,AN,C,m,x8,YS,I0,n 

A fiMAX, JMAX,CM i02,A0S0,DX2,0Y2,KL0SE , 1 1 , J3UP, JSON) 

IF (N.EO.O) RETURN 

MR1TE(4,50)N 

HHITE(A,A0) 

WR?TE(6,«0) (X8(l),I*l,N) 

NRITF(6,70) 

NRITE(6,«0) (Y8(n,I»l»N) 

WRITE (6,40) 

RETURN 






1245 

1246 

1247 
124S 
1245 

1250 

1251 

1252 
125J 

1254 

1255 

1256 

1257 
1256 
1255 
12N0 
1261 
1262 
1261 

1264 

1265 

1266 
1267 
1 ?6B 
;?69 

1270 
12TI 
1272 

1271 

1274 

1275 

1276 

1277 

1278 

1279 

1280 
1 2AJ 
1282 
1281 
1284 
1265 
1286 

1287 

1288 

1289 

1290 

1291 

1292 

1293 

1294 

1295 

1296 

1297 
1296 
1299 
J 300 

1301 

1302 
1301 
1 104 
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40 fOUMAT (//> 

50 FORMAT (TH| Na«ll//) 

40 FORMAT (tSH X8(K), 

70 FORMAT (///J8H YSTK), 

60 FORMAT(1X,8CIO,3) 

END . , 

SUBROUTINE MCHART (R, AK, AN,F ,6» TBR » 0R0»8T,CT»LS, *•» I0» JSXR 
a «tMAX,jHAX,6M]02«A0$0>DX2,0V2>KL0SE,in 

••••••P1.0T8 CMART of local MaCm mumbER, MORITONTAL R0*<S ARE IsCONS 

dimension p;io»n» *K(t)» AN(i)» T(n. Gti), ybr *)» oponi/ sun 
n CT(|), Lsnn m(,j 

RR!TC(6» 170) 

I«l 

DO in Kst * JMAX» J8KF 
iO M(K)nJMA3«l*K 

RRITE(6.200) (M(X) #K«l» J max, JSKP) 
m(JMAX)s100./SQRT(AOSO«CM1O2) 

RR1TF(8, 180) 

20 JajMAK 

If d.EO.n 60 TO SO 
IMRI-1 

|mn«1«2 

IF (I.EQ.2) Imm*2 
F0«F(I)*0X2 

IF (I.CO.lMAX) 60 TO SO 

IP«I*l 

lPPalt2 

IF I^PsI 

SO JP«J*1 

DPJ«L(J)«0v2*(P(1,Jm)»P(1» JP)) 

IF (I.eo.n 60 Tn HO 
IF (I.EQ.IMAX) 60 TO 150 
60 TO bO 

«0 0PIa3.*P(l» 

60 TO 70 

60 0H"P(IP«J)«P(1M.J) 

70 IF (I,f.T.in CO TO 80 

U«CT(naFOaDPI/(l,*AK(I)«AN(J)) 

60 TO 90 

80 Ual.4FD*0Pl» YBP(n*OPJ 
90 V»-8T(!H0PJ 
UUaU«U 
VV«/*V 
00«UUAVV 

AAaAOS0«6M102*0Q 
IF (AA.cT.n.) 60 TO Jon 
wRITE(6,lbO)I,TiUll,VV,00,AF 
RETURN 

100 CONTINUE 


XBJMAXA1*J 




M(K)b100,*8QRT(CO/AA) 

IF U.EQ.2) 
JbJ«J3KP 

GO 

TO 

no 

IF (J,LT,2) 
60 TO SO 

110 continue 

GO 

TO 

no 


isni, 

1306 

1307 
iSnS 

1309 

1310 

1311 

1312 

1313 
1 3t<i 

1315 

1316 

1317 

1318 

1319 
1 3?0 

1321 

1322 

1323 
1 324 
1 325 

1326 

1327 

1328 

1329 
1 330 

1331 

1332 
1 333 
13*4 

1335 

1336 

1337 

1338 
1359 

1340 

1341 

1342 

1343 

1344 

1345 
1 346 
1347 

1346 

1349 

1350 

1351 

1352 

1353 

1354 

1355 

1356 

1357 

1358 

1359 

1360 
1561 

1362 

1363 

1364 
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)«RlTe(6»t90) I.(M(K),K<|,JmAX,JSKP) 1S69 

If (I.EO.IHAK) CO TO 120 lift* 

l«t«l 1S67 

00 TO 20 

120 CONTINUE 1369 

nHITE(6,180) 1370 

00 130 K«l,JN*x»J8*P 1371 

130 M(iO«jMAXtl-K 1372 

HRITE16,200) (N(X)»K*l,JMAx,JSKt») 1373 

RETURN 1370 

1«0 U*0, 1375 

GO TO 90 1376 

150 IF (K'.OSt,EO,0) CO TO «0 1377 

GO TO 190 1376 

C 1379 

C 1360 

160 format (////• negative SPEtO OF SOUNt. OCCURRED In mC^aRT AT» 1361 

1* .*01NT .1*,1«//6M HUB, Ell, 9, 9N VVsHl.u 1362 

2 A Q09*,F1 I ,«»9H AA«,E1 1 ,9////) 1363 

170 F0RM*T(/*1 “ACh no, CH*RT*/1 1369 

ISO FORM*! (/) 1 565 

190 format (19, 9H// ,3119) 1366 

200 FORMAT (6X,3lla) 1367 

END 1366 

SUftROUTl'if SONLlN (P,F,ST,CT,K6,VB,A9,FM,YHP,f jM.AN,r,,M,xs»7S,T0,N 1369 
A ,1m*X.JmAX,GM102,A0SQ,DX2,0Y2,KL0SE,I1,JSuP,JS0K) 1390 

C *191 

c •— -•CAl CULATES XS,VS coordinates OF SONIC LINE — *— 1392 

C 1 593 

OInENSIOn P(TD,n, F(l), STM), CT(l), XOdl, Ybd) AAd)# 1 39« 

1 FJMd), XSd). VSd). YBPd), AN(l), 6d), M(l) 139S 

C 1396 

0(QO)aAOSO-G11C2»OQ 1397 

AM«Ch(RO)sSORT (CfV0(t;U) ) 1396 

C 1599 

00 10 J*1.JMAX luftO 

10 M(J)«0 1901 

NbO 1902 

JaJMAX 1903 

20 1*1 1909 

30 OP1«0, 1905 

GO TO 50 1906 

90 OPla(Pa*l ,J)-Pn«t ,J))aF(1)*0X2 1907 

50 0PJ*(Pn,J-l)»P('.,JAl))*C(J)ADY2 1906 

MR*l,/d,A/.K(l)AAN(J)) 1909 

IF (1,GT,I1) H»el, 1910 

U^CTCDADPl.MR-ORjAVBPd) 1911 

V*«8T{n*0PJ 1912 

OQeUAUAVAV 1913 

OA«0(>)Q) 1919 

IF (OA,CT,P.) CO TO 60 1915 

WRITE ( 6 , 190) I,J,N,P(Ul,J),P(l.l,J),OPl,P(l,J«n.P(l»JAl),CPJ,U.V 1916 
A ,QQ,OA 1917 

RETURN 1916 

60 continue 1919 

FM(I)tAM4CH(f)Q) 1920 

IF (1,E0.1) GO TO no 1921 

Hl*FMn.n«l, 1922 

MtFMin.l, 1923 

HSrHAHI 1929 
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IF (HS.ce.O.) GO TO 100 

1425 


IF (JSON.tO.O) GO TO 70 

1426 


IF (I.EQ.I'IAa.ANO.KLOSE.EQ.O) GU TO 120 

1427 

70 

CONTINUE 

1426 


N«N«t 

k929 


IF (N,LE,396) go to 90 

1950 

SO 

WR1TC(6,200) 

1931 


RETU9N 

1452 

90 

CONTINUE 

1433 


XaX9(n«AN(J)*3Tm 

1934 


XlsXBn-n-AN(J)*ST(I-l) 

1455 


v*v8n)A*N(j5«cT(n 

1436 


Yl«YB(!-n+AN(J)*CT(l-l ) 

1437 



1938 


X3(N)*X1+H2*(5-X1 ) 

1959 


YS(N)»Yl*H2*(Y-Yn 

1990 


M(J)8M(J)+1 

1991 


1F(N,EQ.39S)G0 to 60 

1942 

too 

IF (I.FQ.IMAX) go to 120 

1993 

no 

IslAl 

1999 


IF (I.EO.INAX) GO TO 150 

1995 


GO TO 90 

1990 

120 

IF (M(J),tO,0) GO TO 190 

1997 


J»J«1 

1448 


IF (J.EO.n GO TO 190 

1449 


IF (J.GT.2) GO TO 20 

14S0 


IF (JSON.iQ.i) fiETURN 

1451 


GO TO 20 

1452 

130 

IF {KlOSF.EQ.n GO TO 30 

1453 


OPla(3,*Pn,J5^U,*P(I-l#J)tP(l*2i Jn*Fa)AOx2 

1 u59 


GO TO 50 

1455 

lao 

IF (JSUP.EO.O) RETURN 

1456 


IBI 

1457 

ISO 

J8JMAX 

1456 

160 

Vs.ST(l)+G(J)*0Y2*(Pa,J-l)«P(nJ + l)) 

1459 


QQaVAV 

1460 


DAsD(OQ) 

1461 


IF (DA.GT.O.) CO TO 170 

1462 


RRITE(6.190) l,J,N,P(l, J«l),P(nJfl),V,QO,OA 

1463 


return 

1469 

170 

FJH(J)sAmaCH(00) 

1965 


IF (J.EQ.JmaX) go to 180 

1966 


Hi9Fj«(vt+n«i, 

1967 


HcFJM(J).l • 

1966 


HS«H*H1 

1969 


IF (H3,CE,0,) GO TO 160 

1970 


IF (N,GC,390) GO TO 80 

1971 


NcnaI 

1972 


XcXB(I)-AN(J)*8T(n 

1973 


xi«xBm«ANCJtn*sT(i) 

1974 


H2>«*H1/(m»HJ ) 

1975 


X8(N)sXHH2*(X«xn 

1976 


YS(N)*0, 

1977 


IF (KLO8E.EQ.0) RETURN 

1978 


IF (I.EO.IMAX) GO TO 160 

1979 


IBIMAX 

1980 


CO TO 150 

1961 

160 

J8 J« 1 

1982 


IF (J.GT.n GO TO 160 

1983 


IF U.EO.lMAX) RETURN 

1989 


40 
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T 


T 


T 



oooo no r» ri o o ri o 



miMAX 1485 

00 TO 150 1888 

C 1«»T 

190 KORHAT (/** NeCATlVE SUUARC OF SOUND SPEED CAlCUUATEO IN SUBRO* 1M86 

1*UT1NE SONLlN*//lX,5n,lOEl2,tt/) 1989 

200 FORMAT (/>* NO. Of SONIC PTS. EXCEEDS 396, SONIC PT, CALCUl** 1U90 

t«TI0NS terminated,*) 1991 

END 1992 

0VERL*T(JERRV,1,6) 1993 

PROGRAM ONES 199fl 

COMMON P(81,82) 1995 

common /BL0K3/IMAX,JMAX,0UM(16),KL0SE 1996 

1997 

••.•••halves MFSH SIZE IN BOTH DIRECTIONS AND USES 9TH»0R0ER 1996 

INTERPOLATION TO DISTRIBUTE POTENTIAL OVER new meSM— — • 1999 

1500 

1501 

••••••renumber I*INDEX such that I»0O0«NN0wN P, I^EVENcUNKnOwN P«*. 1502 

1P»1MAX*1 1503 

Ma2«IMAX*l 1504 

00 10 J»1.JMAX 1505 

DO 10 Xai.lMAX 1506 

H1bM«2*k 1507 

M2S1P-K 1506 

10 P(M1, J)aP(M2,j) 1509 

IHAXa2*lHAX>l 1510 

1511 

••••renumber J«INDEX similarly— —o 1512 

lMAXsM.2 15*3 

JPrJHAX*] 1519 

Na2*JHAX*l 1515 

DO 20 191,IMAX,2 1516 

DO 20 K»1,JMAX 1517 

NlfN»2*K 1518 

N2«JP»K 1519 

20 P(I,N1)»P(I,N2) 1520 

JMAXbN«2 1521 

MsIMAX*! 1522 

NeJMAX*! 1523 

1529 

••••••NOW FILL IN 000 J*R0wS#BUT TREAT 1*2 AND laIKAX«l FIRST TO 1525 

ACCOUNT FOR symmetry OR END CONOIT ION.»..— 1526 

1527 

DO 30 J«l,JMAx,2 1528 

30 P(2»J)»,5625*P(1.J)*,5*P(3,J)®.0625*P(5,J) 1529 

IF (KLOSE.EQ.l) CO TO SO 1530 

C 1531 

C ••••••lalMAX IS NOT A symmetry AXIS, SO USE NONCENTRAL INTERP,— 1532 

DO 90 J«1,JMAX,2 1533 

90 P(M,J)*,’^-'5*(P(M*1,J)»P(M*5,J))*,9375*P (('•!, J)*,0625*P(m,5,J) 1539 

CO TO 70 1535 

C 1536 

C ••••••laiMAX IS A symmetry AXIS— •• 1537 

50 DO 60 J«1,Jmax,2 1538 

60 P{M,J)B,5625*PfMtl,J)*,5*P(M«l,J).,0625*P{M«3,J) 1539 

70 m«m*2 1590 

DO 80 J*1,JMAX,2 1591 

DO 80 I«9,M,2 1542 

SO P(I,J)«,5625*{P(I*1,J>*P(I-1,J))-.062S*(P(I*5,J)*P(I-3,J)) 1593 

C 1599 




e •••.••NOW *LL I»IN0ICES are known on all 000 J, KILL IN *LL even J 15«5 


C AFTER FIRST TREATING J«2 AND JMAX»1 BY NONCFNTRAL IntERP, t5«6 

<* J5«7 

00 90 i*t»iMAx isae 

90 P(1.2)«,Jl?5*(P(I,l)-P(l»5n*,9375wP(l,3)^,0625wP(I,7) 1 SP 9 

00 100 lalilHAX ISSO 

100 P(lfN)a,3125*(P(I«N«l).P(l,N.3)}«,9375*P(l,N-l}«,062SwP(I,N.S) 1S51 

Nbn.2 J552 

00 110 X*1,IMAX 1553 

00 l\0 J«a,N,2 155(4 

no P(if j)s,s625*(P(i»j*mP(i,j-n)-,06?5*(Pn,jt5)*p(i,j-5)) 1555 

return 1556 

ENO 1557 

OVEHLAY(JERRY,2,0) 1558 

program TwOO 1559 

dimension XB(200), YB(200), CP(200)i 0ESC(6) 1560 

dimension xs(aoo), Ystaooj i56i 

dimension AN<ino),3T(200),CT(200), 01(200), 02(200) 1562 

call PSEUOn 1583 

call LEROY 158(1 

rewind « 1565 

10 READ («) DESC 1j66 

IF(EOF(<l)) 20,50 1567 

20 CALL CALPLT (0,0,999) 1568 

RETURN 1569 

30 CONTINUE 1570 

read (U) Ihax, jmax,IT,KLOSE,NSL 1571 

READ («)CP5TAR, AMI NF,r>PH,xREF,OXlOxn,DNDYO,OF3 1572 

IF(KLQSE.EO,nGO TO 51 ^575 

REAO(«)CXM,XN,XIM,OXIOXM ,57(4 

31 continue 1575 

READ («) (An(J),Js2, JRAX) 1576 

READ («)fST(I),Ial,lMAX) 1577 

READ (a)(CT(l),l=l,lMAX) 1578 

READ (R) (XB(I),Isl,IM4X) 1579 

READ (0) (YB(I),I=1,1HAX) I5fi0 

READ (0) (CP(n,Isl,lMAX) 1581 

IF (NSLiEQ.O) GO TO <10 1582 

READ («) tXS(l),Isl,NSL) 1583 

READ (tt) (YS(n,I*l,NSL) 158« 

flO CONTINUE 1585 

call GR10(IMAX,Jmax,XB,Y6,ST,CT,AN,01,02,XREF,klOSE) 1586 

CALL PLOT (IMAX,Jmax,KLOSE,XB,YB,CP,OESC,IT,AMInF,CPSTAR,0Pm,xR 1587 

lEF,0XIDXn,DNDY0,CXM,XM,XlM,DXI0XM,nF5) 1 588 

IF (NSL.EO.O) CO T'' 10 1589 

CALL SONPLT (XH, VB,XS,YS,NSL,INAX,XREF,KL08E) 1590 

GO TO 10 1591 

ENO 1592 

subroutine GRI0(IMAX, JMAX,XB, YB,ST,CT, AN,D1 ,02,XPEF,KLOSE) 15P3 

DIMENSION X9(l), YB(l),AN(l), OKI), 02(1 ), 37 ( 1 ), CT(1) 1599 

0XRB5./XHEF 1595 

XSHIFT*3, 1596 

YSHIFTa2,5 I597 

IM«Imax* 1 +klO$E 1598 

00 2 J>'2,<Jhax 1599 

00 1 Ial,lM 1600 

0KI)«(XB(I)-AN(J)*3T(I))*0XR ♦ XSKIFT 1601 

02(l)a(YB(nAAN(J)*CT(n)w(lXP ♦ YSHIFT 1602 

1 CONTINUE 1601 

CALL 0RAW(01,D2,IM) 1609 


42 





i CONTINUE 
00 1 1*1*1N 

D2(n*»V6(n*0XR»YSHIFT »60T 

S CONTINUE ><»0® 

CALL ORArtCOl ,02»IM) 1609 

00 5 lsl,lM ><»10 

00 « J«2,JMAX Ifcll 

01 (J»S)»(XB(n«AN(J)*STn))*OXR ♦ XSHIFT 

02(J-J)«(YB(n4*N(J)*CT(l))*OXR ♦ YSHIFT 1613 

H CONTINUE 

CALL 0»AH(01,D2,JNAX.n *‘‘^5 

S CONTINUE 

CALL NFRA»<E 1617 

RETURN >616 

END >«•>’ 

SUBROUTINE PLOT ( T « AX , JM AX , KLOSE , XB, YB , CP . PESC , I T , Af< I nE , CPST4R # 1620 

1 OPN*XREE*OXIDXO,DNDYO,CXh,xm,XIH,OXIOXM,OE3) 1621 

OlHENSION T(30)»LBLE(6 )»nanE(13) 1622 

OINENSION XB(1), YB(1), CP(l)i DESC(l) 1623 

DATA NREAO/0/ '6?a 

DATA NAME/2hms,7H, 1>^AXs,7H, JHAXs.SH, ITs,fcH, DPMs,7h0x10X0s, 1625 

1 6H, 0N0Y0*,5H CXHa.SH, XMc,6H, XlMs,RH, OXlOXMa,6H CDs 1626 

A ,6K» Qt-5*./ 1627 

NREA0SNREA041 1628 

IF (NREAO.GT.n 60 TO 10 162R 

CALL JPARAMSCn 1630 

YP6»7, 1631 

YDVSO, 1632 

YTlCB.l, 1633 

10 CONTINUE 163(i 

CALL CALPLT(2.0,2,S,«3) 1655 

IMsIMAX-ltKLOSE 1636 

XB(lM+n=0, 1637 

YBnM4l)sX0(lN4n 1638 

XB(!H42)s,2 1636 

Y6(IM4?)sX8(IM42) 16U0 

YMAXtO, 16AI 

OXRPl./XREF 16«2 

DO 20 Ts 1,IM 16a3 

X0n)sXB(n*OXR 16Ua 

YB(l)aY8(n»DXR 16i(5 

IF (YBCn-YNAX.LE.O,) GO TO 20 16^6 

YMAXBYBtl) 16«7 

20 CONTINUE 16UB 

NBOOoO 1609 

IF (YMAX*5,,GE,l,34l,E«06) NBOOa2 1650 

IF (YMAXA5,,6E,2,S4l,r-06) NRODal 1651 

CP(lH4l)al,5 1652 

CP(IH42)a«,5 1653 

0L«CP{IH+1) 165« 

TLaBL4VPr,Arp(lM42) 1655 

PY0a(8L"CPSTAR)/CP(lM+2) 1656 

PYusaPYO 1657 

DO ao 131,1M 1659 

ot(cp(n-cpnH4n)/cp(iM42)42,5 i6S9 

ir (O.LE.IO.) GO TO 30 1660 

CP(n3CP(lM4l)47,S*CP(lM42) 1661 

CO TO ao 1662 

3'3 IF (D,GE,0,) GO TO <'0 1663 

CP(i)aCP(lH4l).2,5*CP(lM*21 1664 


\2 




\ 


1 


t 


40 CONTINUr 

NSVMI22 

NtINE«t 

CALL AXES (-,5,0,,40,,YP6,CP(lM*n*CP(lMf2),YTlC,VOV,2HCP,,20, 

1«2) 

IF (CPST*R,LT,TL,OR,CPSTA«,GT,flL) 00 TO 50 

e 

C DRAH line for cpstar 

c 

CALL CALPLT (-.5,PYU»3) 

CALL CALPLT (,28,PVU,2) 

C 

C PLOT CP 

C 

So CALL LINPLT (XB#CP,IN,j,nlINE#n5YM,i,0) 

IF (NftOO.EO.n 60 TO 70 

c 

c PLOT BODY 

C 

CALL LINPLT (XB«Y6.IN,1»0,0«1»0) 

IF (NBOO.EO.2) go TO 70 
DO 60 

60 YB(ns-Y6(I) 

CALL LINPLT (xe,YB»lM,l»0, 0,1,0) 

70 CONTINUE 
C 

C ADD labels 

C 

•CALL N0TATE(-.5,-t ,ia,orSC.O.,60) 

.ENCOOt(50,aO,LBtE)-i‘^Em, AMTNF,NAME(2.i,lN4)(,NAME(3),J‘'AX,NA«t («), 

in.NAMECS) .DP^' 

call N0TUE(-,5v»ltfc4, ,l«,L0l-C,0,,S0) 

ENCODE(3<»,90,LBLE)NAN£(6),OXIOXO,NAME(7),PNOYO,fJAHF(15) ,QF3 
CALL N0TiTE(-,5,-t.e9, .10,LBLE,0.,39) 

1F(KLOSE,£0.0)GO TO 72 
0Y«>2.Mi 

71 ENCODE ( 17, 100, LBLE)T(U,T (23) 

CALL NOTaTE(-,S,OY, ,14,L8LE»0,, 17) 

GO TO 73 

72 ENCOOE(5fl,llO,LBLE)NAME(6),CXM,NAME(9),XH,NANiE(10),XlM,NAM£(in, 

lOXIOXM 

CALL NOTaTF(-,5,-2.1«, ,1«,LBLE,0,,SO) 

DY*«2,39 
GO TO 71 

73 CONTINUE 
CALL NFRAME 
RETURN 

C 

60 F0RMAT(A2,F5,3,A7,I3,AT,13,A5,ia,A6,E6,2) 

90 F0RNAT(A7,F5.2, A6,E8,2,A6,F5.2) 

100 FORMAT(A7,A10) 

no F0RMAT(A5,E8,2, A5,E6,2,A6,E6,2,A9,F5,2) 

END 

SUBROUTINE 80NPLT (X0, Y6, XS, V$ ,NSL, l^AX, XREF , KLOSE ) 

C 

c •••••-scales and plots body and SONIC LINES — — 

C 

dimension X0(2OO), YB(200), XS(«00), YS(ROO) 

C 

iMalMAX-nKLOSE 


1665 

1666 
1667 
1666 

1669 

1670 

1671 

1672 

1673 
16741 

1675 

1676 

1677 

1678 

1679 

1680 
1681 
16S2 
1685 

16841 

1685 

1686 

1687 

1688 

1689 

1690 

1691 

1692 

1693 
169<1 

1695 

1696 

1697 
1696 

1699 

1700 

1701 

1702 

1703 

1704 

1705 

1706 

1707 

1708 

1709 

1710 

1711 

1712 

1713 

1714 

1715 

1716 

1717 

me 

1719 

1720 

1721 

1722 

1723 

1724 
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T 




T 


T 



o o o r» o o 


r 


I 


DXRM./XREP 
DO 10 I«1»NSL 
XS(n«XS(n*OXR 
10 YS(n«YS(n«oxR 

CALL 09CALE (X9,NSL«X8hax,XSmin) 
call DSCALE (Y8*N8L»Y?max, YSMIN) 
XMAXaXSHAX 

IF (Xfi(lM) jGT.XMAX) XMAXbXB(IH) 
XHir^aO. 

IF CXSHlN.LT.O.) XmInbxSmin 
DO 20 I«liIM 
20 Y8ma488(YP(I)) 

CALL DSCALE (YB,IM,YBM*X, YBMIN) 

OXaXHAX«XMlN 

OYaYSMAX 

Lai 

DXRai, 

30 IF (DX*DXR,LC,2,u>i ,E,06) GO TO fcO 


following card gives further size rfduction if you HE'^OvE COhmENT 

GO TO (00,50), L 


IF(L,E0,2) GO TO 60 
00 0XR3,S*0XR 
La2 

GO TO 30 
50 DXRa.OAOXR 
L = S 

60 IF (OYwDXR.LE. 1 ,5*1 ,E«08) GO TO 00 
GO TO (70,80,90), L 
70 PXR3,5*DXR 
L"2 

GO TO 60 
60 DXRa,a«DXR 
L«3 

GO TO 60 
90 Kao 

00 no lal.NSL 
IF (L.EO.I) go TO 100 
X8(I)=XS(I)*OXR 
Y9(naY3(I)*0XR 

100 IF (YS(I)»1 ,5.GT, I ,E-08) GO TO 110 
KaK*l 


XS(K;aX9(l) 
YS(K)cyS(I) 
no COMTINUF 

IF (L.EQ.l) GO TO 130 
DO 120 1 = 1, 1-^ 
XP(I)sPXR*XH(I) 

120 YB(1 )BOXR*YH(I) 


XHINsXmin.dXH 
YBMAxavtivA X*nxR 


XMAXaXWAX*DX9 
130 XBdM+nsO, 

Y8f IWYnsXRdMAl) 
XS(K»naX6dH*n 
YS(KYnaxBdH*l) 
X0(IM4;)8,2 




1725 

1726 

1727 

1728 

1729 
1750 

1731 

1732 
1753 
175« 

1735 

1736 

1737 

1738 

1739 

1700 

1701 

1702 

1703 
1700 
1705 
170b 
1 /07 

1708 

1709 

1750 

1751 

1752 
1755 
1 7b0 

1755 

1756 

1757 
1756 

1759 

1760 

1761 

1762 

1763 
1760 

1765 

1766 

1767 

1768 

1769 

1770 

1771 

1772 

1773 
1770 

1775 

1776 

1777 

1778 

1779 

1760 

1761 

1762 

1763 
1760 
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— ^ 


T 


I 



1 


I 


I 


ve(iM42)ixenM«2) 

X8(K«2)axe(lM«2) 

V8(K*2)*XB(|M^2) 

8KIPsS,*ARS(XHIN)«2, 

NSKlPiSKIP 

8KIPSN3KIP 

CALL CALPLT(SKlP,2,5#-5) 

CALL LINPLT (XB»VB*IM,l,0#0#l»ft) 

IF (YBMAX*S,.CT,2,Stl.E«08) GO TO ISO 
DO lOO lal.lM 
i«0. Y8(l)«-YB(n 

CALL LINPLT (XB,YB,IM,1»0,0»1,0) 

ISO CALL LINPLT (XS# YS,K, 1,-1 ,22, I .0) 

CALL NFRAMf 

RETURN 

END 

SUBROUTINE OSCALE (X,N,XNAX,XNIN) 

C 

c — — CO'^PUTES MAX AND min OF X-ARRAY- 

C 

DIMENSION X(i) 

C 

XMlNSX(l) 

XMAXMXMlN 
DO 20 181, N 

IF (Xdl.CE.XMIN) 60 TO 10 
XMINaX(l) 

CO TO 20 

10 IF (xm ,LF.,YM4X) go to 20 
XHAXaX(I) 

20 CONTINUE 
RETURN 
END 


17B5 

17B6 

17R7 

1786 

178R 

1700 

1791 

1792 
179J 
179« 

1795 

1796 

1797 
179B 

1799 

1800 
1601 
1802 
1803 
160(1 

1605 

1606 

1607 

1608 

1609 

1610 
1611 
1612 
161 J 
1 61 a 

1615 

1616 
1617 
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APPENDIX D 
SAMPLE CASES 

The input for the sample cases is listed below. The output for these 
cases is on the following pages. Note that the sample cases are for only 
25 cycles on the crude grid. In actual usage there might be more cycles 
and some grid refinements. The plotted output is shown in figures 1 to d. 

10«l CLllPSOIO 


5S 

0,000000 

0,012311 

0.038155 

0,072396 

0,113710 

0,161530 

0,21521 

0,27«21 

0,33799 

0,90601 

0,97777 

0,55279 

0,63C«l 

0,7102* 

0,79181 

0,8;-Sr' 

0,9S709 

1,09190 

1.1296 

1,2069 

1,2876 

1,3668 

1,9935 

1,51 73 

1.S077 

1,6593 

1.7169 

1.7737 

1.82S5 

1,8715 

l.PllO 

1.993e 

1 ,96AS 
,000000 

1,967689 

,015693 

2 f 

,027360 

,037399 

,096312 

,059995 

,061977 

,068792 

,07«9«9 

,080997 

,085280 

,069690 

.092OJ6 

.095711 

,097609 

,099210 

,099911 

,099019 

,099221 

,097*37 

,095770 

,093026 

,089627 

,0*55*0 

,000906 

,075626 

,069768 

,063360 

,056959 

,099097 

•09l2tt1 

,033100 

,0?a7tS ,0SS6itl 

999, 999, 

?1 21 2S e 

!•« 1.0 

.5 1,3 

l,a ,9«S 

tO-t tUTPOOlO -<ITm 
19 

0,0000 0,012311 

0. 

.1 2. 

' 1 1 
o« 

.089 

?0*PE0Ct6l S716C 
0,03*155 0,072396 

0,113710 

0, 16153P 

0.21521 

0,27**21 

0,33799 

0,90601 

0,67777 

0,5S27u 

0,63?9l 

0.71027 

0,79181 

0.67C52 

0.9S7S9 

1,09190 

1.1296 

1,2069 

1,2876 

1,3668 

1,9935 

1,5173 

1,5077 

1.6593 

1.7169 

1.7737 

1,8255 

1,8715 

1.9110 

1,9936 

1,9688 

,000000 

1,9797950 

,015693 

,027360 

,0373116 

,096312 

,059u95 

,061977 

,066792 

,079999 

,080997 

,085260 

, 0 ^ 91.90 

,092CJO 

.09571 1 

,097*09 

,099210 

,099911 

,099919 

,099221 

,097337 

,095770 

,093029 

,089627 

,085560 

,080906 

,075626 

,069768 

,063360 

,056939 

,099097 

.P9J291 

,033100 

,020765 

9999, 

21 21 

1.8 

.5 

.020000 
•,98989795 
25 0 

1.0 

1.3 

.1 

0 

0. 

,C89 

2. 

1 

1,9797959 

• >3 

2 . 




SPHiPE/iS-ofc cnKF/r»tl*ii>f«/i‘i-rf 6 rtiRF 

4*0 

,0 7,005F-05 1,751F«C2 S.SKuF-ra 

7,lO»-F»ri l.insf^oo 1,5'?€400 !,5?6f400 

2,l?3Ui.O 2,7jeC4P0 5,ce^»*4f'0 >,7S5F4«0 •,072F#"‘ft 

а, «0frF4nD ii,7^«i4nCl V,fiS5FOC 

б, ^«aF♦nO 6,9lO»4Pf 7,li^t#oe T,\S?tfPC 7,»,««^4rp 7,BS#»t4ftO 0,IS6F4*>0 ft,sn^4''0 

,0 U512F-PI \,bbui*r^\ ?,<i67f«ftl 

<i,9n''F-A| S,29ttf*CI *>,7-lt-ni 

6,77iF*0| l.OPOtano I.P^-OFa'^O 140^0^4^0 

j,o?er40o i, 202^400 i,?<‘<.F40o j,3M»4f.o i,«5#*moo i,so9f4flo 
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54 

1,979796 

,C2C0CC 


OYDk^i* 99R</,00e0 
DYDXT» -.weRV 
tHAi» ,1000 

MRtf* 2,0000 


52 


T- 


T 


T 


r 



|MA«l i\ 

JHAXA 21 

Mlt* 2% 

NMAUr«0 
KIOSCAO 
NPlUYal 
RFiaifMOO 
COViRUa ,lOt «01 
OFJ» 0. 

ONOYOa ,S0Ot«O0 
ALFal ,30 

DXiDAOa , 6 A 0(»01 
XNa 

CXHa • 7 SOUOU 
OXlUXha , 200 t 401 
CAHBl.aO 
AMiNFa ,9950 


1 

s 

t 

7 

THIT 


AK 

f 

X 

0. 

0* 

0 • 

,90001402 

,90001402 

, 6967fc 4 v2 

, 1 I90t 402 

2 


,156lt-02 

,56791-02 

.59211402 

,59211402 

,656«t 402 

,**S67t4y 1 


,2609L*01 

, l699t-0l 

, 16571-01 

,276U 402 

,27611402 

,1161f402 

, 165U 4ul 

4 

,7056L«01 

,5B4 ifOl 

,35blt-0l 

,l555t402 

, J 5551 40 2 

•25261 401 

,r297t *oC 

•» 

, ia9oe^oo 

,15511400 

,50201-01 

,97651401 

. 97651401 

, 7 j9 5t 4 cG 

,51 22t 400 

(> 

,2i.feSt4 00 

,25l6t400 

,66521 -01 

,6‘.3 if 401 

,69311 401 

,5361140^ 

, i 5 5 L ♦ 0 C 

7 

,923St4C)0 

,90621400 

,«0Llt-01 

,91991401 

,91991 401 

, 1 66 ff 4t f 

,2'“«Gl4wC 

6 

,bU st400 

,6U09t 400 

,9lb7t -01 

,299t,t 4Ul 

, 29 '.tf 4V« 1 

, ‘ ’ J t 4 0 f 

, 2 - 1 1 k 4 w C 

9 

,»i»*t400 

,SlS7l4CiC 

-CM 

< n K / l « j 

, 10571 401 

, i U^>U 41 0 

, 2. t f ^ 

10 

, 1 0febl401 

, iOS0l40l 

,996»‘.i-0 1 

oi 

• ,269 il 40( 

, 1 f^OcL 4 t 0 

, f 1 56t ♦ ' ' 

11 

,1? it*Ol 

,127 ft 401 

, 96C9C-0 1 

0, 

• ♦ 1651 1 401 

, 1 1 52t 41 0 

, 7 wt 4 W i 

12 

, ISCOt 401 

• 1 U69t401 

,67501-01 

0, 

^,31671 401 

. 1 5 0 it 4 c 

,c^c2t4v'G 

15 

, 167iik 401 

,16561401 

, 755«t-01 

0, 

-,99911 40 1 

,2l/6t4ci' 

. 325M4‘ 

lA 

,10071401 

,17911401 

,61201-01 

0, 

-,73721 401 

,9Jl6t40C 

, - SS wt 4 V C 

IS 

, 1906tUi 

,169 01 *; 1 

,9562t-0l 

0. 

- , 11 05l 412 

, 1 6 0 L 4 0 1 

,5*^t2t4u{ 

d 

, I99bt 4ul 

,19^01401 

,20^01-01 

0, 

•,2610f 402 

, 125‘^L4o? 

,5 'Wf't4wt 

»r 

, 2 1 2 U 4 y 1 

,2U'5t4. 1 


0, 

0, 

0, 

, 321 C f 4 V ( 

IS 

,23296401 

,25131401 

,200C1-01 

0# 

0. 

0, 

• i-^at^wr 

19 

,27961 401 

,27 .01401 

, 2 0 0 0 1 • 0 1 

0. 

0. 

0, 

, t y r t - w 1 

20 

, J996t4ul 

, 59r0t 401 

,2>0Ct-Cl 

0. 

0, 

0* 

, 2 L- 0 C t - t 1 

2: 

, I000t4>! 

,10001451 

,20;'^1-01 

0, 

0. 

0, 

, 1 f vGL-c 9 
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T 


T 


T 





J 


I 


NOHHAt COUfiO, STRETCH FOR ALFb 1,}Q0 


J 

AN 

c 

6H 

1 

0. 

0 • 

,79206-03 

2 

,233CE^02 

,158«£»02 

,47366-02 

3 

,6979E*01 

,7693E«02 

,10096-01 

a 

,5006E*01 

,2030E-01 

,30056-01 

5 

,32<iie«0l 

,39fllt-0l 

,53576-01 

6 

•2274E401 

,6732t-01 

,85«9E-0l 

7 

,1674£401 

.1037E400 

,12666400 

e 

,1272E401 

, 1496C*00 

,17786400 

9 

,9873E*00 

,2060L400 

,23986400 

10 

,7765E+00 

.2736E400 

,31346400 

u 

,6156t+00 

,35326400 

,39936400 

12 

,d695E+00 

,dU55E400 

,49856400 

13 

,36ebt*00 

,55156400 

,61186400 

\9 

,3:64E*00 

,07206400 

,73996400 

15 

.2365E400 

,60786400 

,08396400 

16 

,iei7fcf00 

,9600t400 

,10456401 

17 

,l357fc400 

,11296401 

,12236401 

16 

.9264E-01 

,13176401 

, 14206-' 1 

19 

,573«£-0l 

,15246401 

,16386401 

20 

.2672E-01 

,17516401 

.18766401 

21 

0. 

,20006401 

,21366401 


CPU SECONDS FOR BODY GEOMETRY CU*^PU T A T I QNSs ,1J2 


IT 

DPMAX 

10 

JO 

RMAX 

IK 

JR 1SU8 ISUP 

Ri wTG 

RFl 

0F3 

•1S 

StC/CTC 

1 

,9146-02 

4 

21 

,939640<f 

1 

21 

0 

,627f 400 

1 , U 0 

0,000 

0 

,209 

2 

,5796-02 

15 

21 

,7036403 

1 


0 

,2026401 

1 ,400 

0, flUO 

55 

.227 

3 

,5196-02 

2 

21 

,5726403 

1 

21 

0 

, 1666401 

1 .4"0 

0 , ooo 

101 

.237 

4 

,5776-02 

1 

21 

,9506403 

1 

21 

0 

,2766401 

1,400 

0 , n 0 ■) 

1 13 

.239 

5 

,4176-02 

1 

21 

,6786403 

1 

21 

0 

,1996401 

1,400 

0,000 

1 14 

.242 

6 

,2656-02 

1 

20 

,2566^03 

1 

21 

0 

,8706400 

1 ,40U 

0,000 

115 

.230 

7 

,1356-02 

1 

20 

,6066402 

1 

21 

0 

,3216400 

I. 4 OO 

0,000 

1 14 

.2ub 

8 

,1386-02 

7 

21 

,79dE402 

1 

21 

0 

,3346400 

1,400 

O.OoO 

113 

,244 

9 

,1206-02 

7 

21 

,1436403 

1 

21 

0 

,4746400 

1 .600 

0,000 

1 13 

.247 

10 

, 1046-C2 

7 

21 

,1546403 

1 

21 

0 

,4896400 

1,400 

0,000 

113 

. 24 7 

1 1 

,6996-05 

7 

21 

,1226403 

1 

2) 

0 

,3976400 

1,600 

0, OvO 

113 

.24 1 

12 

,7786-03 

7 

21 

,6606402 

1 

21 

0 

,300640C 

1,400 

0,000 

112 

.245 

13 

,6756-03 

7 

2! 

,7086402 

1 

21 

0 

.2476400 

1,400 

0,000 

1 l‘3 

, 24 u 

14 

,6016-03 

b 

21 

,6596402 

1 

21 

0 

,2266400 

1,400 

0,000 

114 

.240 

15 

,5396-05 

8 

21 

,6256402 

1 

21 

0 

.2106400 

1,400 


112 

.246 

16 

,4616-03 

8 

21 

,564640^ 

1 

21 

0 

,1696400 

1,400 

0,000 

113 

,247 

17 

,4306-03 

6 

21 

,4916402 

1 

21 

0 

, 1666400 

1,400 

0,000 

113 

,242 

16 

,383t» 3 

8 

21 

,4306402 

1 

21 

0 

,1466400 

1,400 

0,000 

113 

• 242 

19 

,3416-03 

6 

21 

,3846402 

1 

21 

0 

,1316400 

1 . 4 OU 

0,000 

113 

,243 

20 

,3046-03 

6 

21 

,3496402 

1 

21 

0 

, 1186400 

1 ,600 

0,^00 

114 

.247 

21 

,2716-03 

B 

21 

,3166402 

1 

21 

0 

,1076400 

1,400 

0, OuO 

114 

,246 

22 

,2426-03 

9 

21 

,2856402 

1 

21 1 

0 

,9596-01 

1 ,<iOO 

0,000 

U4 

,26(^ 

23 

,2216-03 

9 

21 

,2576402 

1 

21 ] 

0 

,8626-01 

1,400 

0,000 

114 

,240 

24 

,2006-03 

9 

21 

,2326402 

1 

21 : 

0 

,7776-01 

I. 4 OC 

0, OuO 

114 

,241 

25 

,1826-03 

9 

21 

,2106402 

1 

21 ! 


,7046-01 

1 ,400 

0,000 

ll« 

,244 

-•-DIO 

NOT converge 

IN 

25 cycles, 

m m m m 

RHAXS 

,216402, 

. COYHb , 

256.02 




CPU S6rONOS# 6 

10 

FOR 

25 ITtRATIDN'S, 

KhALFrI 
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f 


2 

IB 

X6 

YB 

CP 

H 

\ 

0,000 

0,000 

0,000 

1,27261 

0,00000 

2 

,006 

,002 

,00k 

1,16736 

.21399 

3 

,026 

.017 

• OIB 

,56402 

•kS98« 

<1 

,071 

,058 

,039 

,20926 

,67617 

S 

,t99 

.135 

,050 

,03982 

.«71«0 

6 

,267 

.252 

,066 

•,C2399 

1,00902 

r 

,020 

• «oe 

,011 

•,05237 

1,02606 

s 

,616 

,600 

,092 

•,06930 

1,03675 

9 

,835 

,619 

,096 

•,07896 

1,09265 

10 

1,066 

1,050 

,100 

-,08511 

1,09691 

It 

1,293 

1.277 

,096 

•r09U4 

1,05011 

t2 

1,500 . 

1,909 

,067 

•,09826 

1,05999 

13 

l,67tt 

1,656 

,075 

•,06619 

1,03963 

ttt 

1.607 

1.791 

,061 

,00611 

,99136 

IS 

1,906 

1,690 

,096 

,12920 

.92216 

u 

1,996 

1,980 

,020 

,38759 

,77228 

17 

2,121 

2,105 

,020 

,11500 

,92798 

te 

2,329 

2,313 

,020 

,02152 

,98221 

19 

2,796 

2,730 

,020 

,00060 

,99219 

20 

2t*A« 

3,996 

3,980 

,020 

,020 

,00006 

0,00000 

,99972 

,99500 


drag COEFFICUM UY TRaPEZOIOAU INTEGRATIONS ,03736 

DRAG COtFFlCUNT BY SJWHSON INTEGRATIONS •,03<IO6 
Plot of cp at uneouau iNCKeMtMs 


I 

XB 

YB 

CP 





1 

0,000 

P,000 

1,2726 

♦ 



• 

2 

,002 

, 006 

1.1679 

♦ 



* 

3 

.017 

.016 

,5699 


♦ 


• 

9 

,058 

,039 

,2093 




♦ t 

5 

.135 

,050 

,0398 





6 

.252 

, 066 

•,0239 





7 

,908 

• 061 

•,0529 




t ♦ 

6 

,600 

• 092 

• •0693 




* ♦ 

9 

• 819 

,098 

•,0790 




• ♦ 

10 

1,050 

,100 

• ,0651 




• ♦ 

11 

1,277 

,096 

•,0911 




* ♦ 

12 

1,469 

,067 

-,0983 




• ♦ 

13 

1.658 

,075 

• ,0661 




• ♦ 

19 

1.791 

.061 

,0061 




♦ 

15 

1,890 

,046 

,1292 




♦ * 

16 

1,980 

,020 

,3675 



♦ 

* 

17 

2,105 

,020 

,1150 




♦ • 

16 

2,313 

,020 

,0215 




♦ 

19 

2,730 

,020 

,0098 





20 

2lt». 

3,960 

,020 
• 020 

,0005 

0.0000 




4 

♦ 




hACH NO 


CHART 



21 

20 

19 

16 

17 

16 

15 

14 

13 

12 

11 

10 

9 

6 

7 

6 

5 

4 

3 

2 

1 

\u 

0 

66 

91 

94 

95 

96 

97 

96 

96 

96 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

in 

21 

65 

90 

V3 

95 

96 

97 

97 

96 

96 

96 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

in 

65 

60 

6b 

69 

92 

93 

95 

9b 

97 

97 

96 

96 

98 

99 

99 

99 

99 

99 

9Q 

99 

99 


67 

67 

69 

90 

91 

92 

93 

94 

95 

96 

96 

97 

96 

98 

96 

99 

99 

99 

99 

99 

99 

5// 

97 

9b 

95 

95 

95 

95 

95 

96 

96 

96 

96 

97 

97 

96 

98 

98 

99 

99 

99 

99 

99 

6// 

100 

100 

99 

99 

99 

96 

96 

9B 

96 

96 

96 

96 

96 

96 

98 

98 

99 

99 

99 

99 

99 

7// 

102 

102 

101 

101 

101 

100 

loo 

100 

100 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

9// 

103 

lo:- 

102 

102 

102 

102 

101 

101 

101 

101 

100 

100 

100 

loo 

99 

99 

99 

99 

99 

99 

99 

9// 

106 

103 

103 

103 

103 

1 u2 

102 

102 

102 

101 

101 

101 

loi 

100 

lor 

loo 

99 

99 

99 

99 

99 

10// 

104 

104 

104 

103 

103 

103 

103 

102 

102 

102 

102 

102 

101 

101 

101 

100 

100 

luO 

99 

99 

99 

n// 

105 

104 

104 

104 

103 

103 

103 

103 

103 

102 

102 

102 

102 

lOl 

lOl 

101 

100 

100 

99 

99 

99 

12// 

lOS 

105 

104 

104 

104 

104 

103 

103 

103 

lOi 

102 

102 

102 

102 

101 

101 

101 

100 

99 

99 

99 

i>// 

103 

102 

10? 

102 

102 

101 

101 

lOl 

101 

101 

lUl 

) 01 

100 

100 

ICO 

IC'O 

100 

100 

99 

99 

99 

i«// 

99 

9ft 

96 

97 

97 

97 

97 

97 

97 

97 

97 

94 

96 

96 

98 

98 

99 

99 

99 

99 

99 

15// 

92 

92 

92 

93 

94 

94 

95 

95 

96 

96 

96 

97 

97 

96 

96 

98 

99 

99 

99 

99 

99 

u// 

76 

64 

67 

69 

91 

92 

93 

94 

95 

95 

96 

97 

97 

97 

96 

98 

99 

99 

99 

99 

99 

17// 

92 

93 

94 

95 

95 

9b 

96 

96 

97 

97 

97 

97 

96 

98 

98 

98 

99 

99 

99 

99 

99 

16// 

96 

96 

96 

98 

98 

96 

96 

96 

98 

96 

98 

96 

98 

96 

98 

99 

99 

99 

99 

99 

99 

19// 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

20// 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

21// 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 

99 


21 

20 

19 

16 

17 

Ifc 

15 

14 

13 

12 

11 

10 

9 

6 

7 

b 

5 

u 

3 

? 

1 


NS Jb 


XS(K), KSl^ ,, 


,224t^00 
,31 UAOO 
,3b7E*00 
,469tfOO 
,759Et00 

, I7bt90l 
, I72t40l 
. I7lt+01 
,l7ffc+0l 
♦ 17U401 

,244E+00 

,32bt+00 

,3P7E+C0 

,502t»0C 

,925E+O0 

,l7St+0l 
, 1 72E + 01 
, 171E + 01 
, 17CE + 01 
, lb9E40l 

,2b9E*00 
, i^OEtOO 
, 41 4E*00 
,549t *00 

, 174E + 01 
i 1721 + 01 
• 17U*01 
,1701401 

,292t400 
,353E4C0 
, 44 1 E 4 OO 
•631E400 

•173E401 
, 17 Ifc401 
, 1701401 
, 17U^01 

tSC<), KSl, 

» t 1 1 » 







,625fc*0l 
,20bE«00 
,4b7fc + 00 
,107tt0l 
.334E+01 

jbuOE^Ol 
, 202f tC.w 
, uSfcEiOO 
, 1 0 0 E ♦ 0 1 
,33lfc401 

.922E-01 
,25bEtOO 
,570EtOO 
, 1 36E+01 
,510EtOl 

,927t»0l 

,250E^OO 

,5bOE900 

,13^Et0l 

,5C6Ef0l 

,1251400 
, 5 1 400 
,b97EtOO 
, 1 7 dE401 

, 124E*00 
,3C6ttO0 
,68bE400 
, l74t4Cl 

, IbiE^OO 
,3S3E-»00 
,8b0E*00 
,237Et01 

, 1601 400 
, 37 oE 4 0 0 
,647Et00 
,23 uE40 1 


CPU SECONDS TO COMPUTE AND PLOT CP AND NCHARTa ,i63 


56 



1 


r 


1 


1 


I 


|PHCRC/l)*0C6 CONE/CYllNDeP/lS»DLG PlARE 


INPUT C00RD1NATE3 
1 X Y 


l 

OiOOOOOO 

0,000000 

2 

•001637 

,040430 

3 

• 007003 

,063390 

• 

• 017510 

,131200 

S 

•035640 

,165400 

4 

•06S090 

,246700 

7 

•110700 

,313600 

4 

•176200 

,362700 

9 

•272700 

,445300 

10 

•396900 

,490000 

u 

•543600 

,529400 

12 

• 710600 

,574100 

13 

,697600 

,624200 

19 

1,105000 

,679700 

IS 

1,332000 

,740600 

14 

1,576000 

,606600 

17 

1,642000 

,677300 

IS 

2,123000 

,952400 

19 

2,421000 

1,000000 

20 

2,736000 

1,000000 

21 

3,065000 

1,000000 

22 

3,396000 

1,000000 

23 

3,735000 

1,000000 

24 

4,072000 

1,000000 

25 

a , *06000 

1.026000 

26 

4,734000 

1 , 1 itooo 

27 

5,053000 

1,202000 

26 

5,361000 

1,264000 

29 

5,655000 

1,363000 

30 

5,9 • aOOO 

1,436000 

31 

6,196000 

1,509000 

32 

6,446000 

1.575000 

33 

6,664000 

1,639000 

34 

6,910000 

1,699000 

35 

7,133000 

1,759000 

36 

7,357000 

1,619000 

37 

7,594000 

1,663000 

36 

7,656000 

1,953000 

39 

6,156000 

2,033000 

40 

6,511000 

2,074000 


DYOXNc 999,0000 
PYOXTb ,20SO 

YNAXs 2,0000 

XPEF* e.sooo 

IHAXc 21 

JNAXa 21 
HIT* 25 
HHALFbO 
KtOStBO 
NPLOTsl 
PFl*l,aOO 
covERGt tioe^oi 
QF5* 0, 

ONOYO* ,266E401 
AtF»l,50 

PXIUXQ* ,200Et01 
XM* t600t40l 
CXH* ,750Ef00 
OXlDXh* ,lS0Et02 
GAH*1,40 
AHINF* «6000 
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1 

1 

1 

i 

8 

! 

1 

t 

I 

1 

i 

1 

Y 

TH6T 

TH6T8 

I 

1 

AK 

1 

1 

i 

p 

1 

0. 

0. 

® f 

,90006402 

,90006402 

,19996401 

,50006400 

2 

,1074t400 

,11486^01 

,10656400 

,77716402 

,77716402 

,20006401 

,40986400 

5 

•2580E400 

,65086-01 

,24676400 

,60466402 

,60466402 

,19946401 

,26956400 

4 

, 49031400 

,22176400 

,41566400 

,35936-02 

.33936402 
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Figure 2. - Frame 2 of plotted output. 
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Figure 3. - Frame 3 of plotted output. 
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Figure 4. - Frame 4 of plotted output. 
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Figure 5. - Frame 5 of plotted output. 
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Figure 6. - Frame 6 of plotted output. 



Figure 7. - Frame 7 of plotted output. 
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Figure 8. - Frame 8 of plotted output. 
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